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COMPUTER PROGRAM FOR POST-FLIGHT EVALUATION OF A LAUNCH VEHICLE 
UPPER STAGE ON-OFF REACTION CONTROL SYSTEM 

SUMMARY 


This report describes a FORTRAN IV coded computer program for post-flight 
evaluation of a launch vehicle upper stage ON-OFF reaction control system. 
Aerodynamic and thrust misalignment disturbances are computed as well as the 
total disturbing moments in pitch, yaw, and roll. Effective thrust 
misalignment angle time histones of the rocket booster motor are calculated. 
Disturbing moments are integrated and used to estimate the required control 
system total impulse. Effective control system specific impulse is computed 
for the boost and coast phases using measured control fuel useage. This 
method has been used for more than fifteen years for analyzing the NASA Scout 
Launch Vehicle second and third stage reaction control system performance. 

The computer program is set up in FORTRAN IV for a CDC CYBER 175 system. 
With slight modification it can be used on other machines having a FORTRAN 
compiler. The program has optional CALCOMP plotting output. With this option 
the program requires 19K words of memoiy and has 786 cards. Running time on a 
CDC CYBER 175 system is less than three (5) seconds for a typical problem. 
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1.0 INTRODUCTION 


Post-flight analysis of a reaction controlled launch vehicle upper stage 
should include computation of the disturbing pitch, yaw, and roll moments, 
control fuel consumption, estimated control fuel specific impulse, as well as 
other performance measures. The method presented herein, provides a useful 
tool to aid in these tasks for an ON-OFF reaction controlled vehicle. It is a 
straightforward deterministic approach based on uncoupled rigid body equations 
of motion. Impulse and control fuel useage is based on integration of the 
absolute value of average duty cycles required to balance disturbing moments 
and optional inputs for additional impulse such as used during limit cycle 
operation. Being straightforward and simple this important evaluation method 
is sometimes overlooked in post-flight evaluation of reaction-control systems. 

The method has been used for more than fifteen years for evaluating the 
ON-OFF hydrogen peroxide reaction control systems on the second and third 
stages of the NASA SCOUT (Scientific Orbital Utility Test) Vehicle. 

Use of this method has led to an accurate data bank on the booster rocket 
motors thrust misalignment and roll torque characteristics and control system 
specific impulse deviations. It has also been useful in helping diagnose 
rocket motor failures. 

Input data required includes dynamic pressure, angles of attack and 
sideslip, angular accelerations or angular rate trace slopes, booster thrust 
and weight time histories, mass properties, control fuel or pressurization 
variables, and vehicle control moment ams. Optional input includes 
incremental rates at control motor firings and/or incremental impulse. 

Methodology and a detailed computer program description is presented 
herein. Details of the computer program including a sample problem and 
detailed input and output descriptions are presented in Section 3* A complete 
FORTRAN listing is presented in Appendix A. 
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2.0 METHODOLOGY 


This section presents the methodology and equations which can be used to 
compute pitch, yaw and roll disturbing moments, effective thrust misalignment, 
total impulse expended by the control system, control fuel consumption and 
effective overall control system specific impusle. Knowledge of vehicle 
angular rate time histories and control motor firings from telemetry data is 
required. Booster thrust, dynamic pressure, angles of attack and sideslip, 
are also required. If the ON-OFF control system is a regulated pressure 
system an option of computing fuel useage is included based on measured 
pressure and temperature of the unregulated pressure supply. 


2.1 Assumptions 


Major assumptions and approximations made in the method are: 

. non- spinning three-axes stabilized vehicle, 

. ON-OFF reaction control motors which are off sufficiently long to 
define angular accelerations about each axis, 

. aerodynamic and mass properties symmetry in pitch and yaw, 

. no gyroscopic cross-coupling terns, 

. aerodynamic coefficients are non-linear with angle of attack and can 
be described by a cubic for normal force coefficient and a linear 
variation of aerodynamic center with absolute value of angle of attack, 

. control impulse expended assumes balance of disturbing moments plus an 
additional impulse supplied as input, 

. control impulse calculations assume independence of pitch, yaw, and 
roll control motors (no mixing of pitch-roll or yaw-roll, etc.), 

. optional calculation of control fuel useage from unregulated pressure 
and temperature assumes nitrogen pressure regulated system without 
venting and non-varying temperature on regulated side, 

. no rocket motor jet damping in pitch, yaw, and roll, 

. no aerodynamic damping 
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2.2 Equations 


2.2.1 Disturbing Moments and Thrust Misalignment 


The angular equations of motion in pitch, yaw and roll baaed on the 
previously mentioned assumptions are (see Figure 1 for sign conventions): 


(2-1) 

lyd = 57.3 2M = 57.3 [m^. + Mr + Mgero] 


(2-2) 

Iy<^= 57.3 2N = 57.3 [nc + + Nae^o] 


(2-3) 

Ix^ = 57.3 2L = 57.3 [tg + L^j 



Control momenta are: 


(2-4) 

Me = -Fcp Ic 

(ft-lbs) 

(2-5) 

Nc = -Fcy Ic 

(ft-lbs) 

(2-6) 

Lc = -2 Fc[^ Rc 

(ft-lbs) 

where the control moment arm is. 


(2-7) 

= (xg - Xgg) cos + Zg sin Xg / 12 

(ft) 


Booster induced momenta which are produced by all sources (i.e., angular, 


offset, swirl, jet damping, etc.) are all lumped into an effective thrust 
misalignment angle in the pitch and yaw planes. 

(2-8) Mr = T Ir / 57.3 (ft-lbs) 

(2-9) Nr = T Cry Ir / 57.3 (ft-lbs) 

where the thrust side force is assumed to act at the nozzle throat, 
(2-10) 1^ = (x, - Xcg) / 12 (ft) 

aerodynamic moments are, 

(2-11) Magr.Q= Cf^ S Q ( x^g - x^p) tana/ 12 tanTj (ft-lbs) 

f^aero=-CN S Q ( Xcg - x^p) tan/3/ 12 tan tj (ft-lbs) 
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total aerodynamic angle (r/) is. 


(2-13) 7 } = tan“^y tan^ a + tan^ /3 


aerodynamic normal force coefficient is a cubic function of ' 77 ', 


(2-14) 



and 

aerodynamic center is. 


(2-15) 

X = X + — P hi 

Xcp ^ da ' ' 

inches 

Total disturbing moments are computed 
when the control motors are off, i.e.. 

from equations 2-1, 2-2, and 2-3 

(2-16) 

Mq = ly e / 57.3 

(ft-lba) 

(2-17) 

Nq = ly ^ /57.3 

(ft-lbs) 

(2-18) 

Ld = /57.3 

(ft-lbs) 

An effective moment due to the rocket 

booster is. 

(2-19) 

Mr = Mq - Magf.0 

(ft-lbs) 

(2-20) 

~ ^aero 

(ft-lbs) 

An effective thrust misalignment can then be computed from equations 
(2-8), ( 2 - 9 ), and (2-l6) through (2-20) for periods when the control motors 
are off. Pitch, yaw and total effective thrust misalignment is, 

(2-21) 

e^p = [^y ^ -57.3 MgeroJ / 3" 

(degrees) 

(2-22) 

^Ty ~ -57.3 Nggy-gj / T If 

(degrees) 

(2-23) 


(degrees) 

Total thrust misalignment angle and orientation angle (\j) is shown in 
Figure 1. 

(2-24) 

Xj. = tan'^ (€,/€,) 

y "p 
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2.2.2 Total Impulse Expended 


Total impulse expended by the control system is the integral of the 
absolute value of thrust-time histones of each control motor. 

(2-25) Ij = E/I F^'jdt ; i= 1,2,3, n motors 

This information is not usually available. Estimation of this from motor 
commands and preflight measured thrust levels is often inaccurate when short 
pulse widths are commanded. During a boost phase the disturbing torq^ues are 
generally large enough for the control system to operate on one aide of the 
deadband resulting in a long period of single direction control motor pulsing 
(Figure 2). When this occurs the average angular impulse expended by the 
pitch yaw and roll control motors balances that produced by the disturbances. 
The linear impulse is, 

(2-26) It=/[|Vc| + |Nd/1c| + |l-D/Rc|]^t + Mj 

where is an additional impulse. 


When the control system deadband is crossed opposite motors fire (Figure 
2) . If this occurs the impulse expended includes the balance of disturbances 
plus the balancing of the opposite motor plus the impulse expended by the 
opposite motor. For this reason the program includes an input table of an 
"additional" impulse (Al) time history which can be estimated by other means 
such as angular rate changes and measured motor pulse widths. 


An option for computing the impulse from supplied changes in pitch, yaw 
and roll rates is available to cover periods of no disturbance such a limit 
cycle operation during a coast period. A typical rate trace is presented in 
Figure 3* The linear impulse required to produce the motion is, 


(2-27) AIt =E 


|lyA0/57.3 Icl + |lyAi|r/57.3 iJ + |ljc^<^>/57 .3Rc| 


With this option the absolute values of incremental pitch, yaw and roll 
rates must be supplied as a time history as described in the input data 
description (paragraph 3*4). 


2.2.3 Control Fuel Useage 


Measured control fuel useage can be input directly to the program to be 
used to compute an effective specific impulse. 


(2-28) I » 


(Ibf-aec/lbjj) 


This IS computed as a running time history during the 
phases. 
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boost and coast 



If a closed pressure regulated liquid monopropellant control system is 
used there is an option to compute the liquid expelled from measured 
unregulated pressure and temperature. The follovfing assumptions are made: 

. total mass of pressurizing gas does not change 
• regulated aide gas pressurant is the same temperature as the 
unregulated aide 

. liquid monopropellant density is constant 

. the compressibility factor of the gas on the regulated side does not 
change 

The euaount of propellant expended is: 

(2-29) Creg Vg ( Pgo ~^g _ K i ^9 

Preg v Cgg Tg^ Cg ) 1 Tg^ 

where , 

weight density of the liquid fuel (ibs/in^) 
compressibility factor of the regulated side gas pressure 
unregulated side pressurizing gas volume (in^) 
regulated side gas pressure (psia) 
unregulated side gas pressure (psia) 
unregulated gas compressibility factor 
unregulated gas temperature (°R) 
initial volume of gas on the regulated side (in^) 

The zero subscript denotes values at the start of operation when a known 
quantity of fuel is in the tanks such as at initial operation. 

The computer program includes a short subroutine (C) which computes the 
compressibility factor for dry diatomic nitrogen baaed on a curve fit. The 
curve fit is accurate for pressures in the range of 400 to 3000 psia and 
temperatures between 50 and 100 degrees Fahrenheit. The curve fit is, 

(2-30) Cg = 0.9977 + 1.657738 X lQ'^(Pg-140Q) + 1.264881 X 10"®(Pg-14Q0)^ 

+ 4.5833 X 10'^(Tg - 522)2 

Control fuel remaining is the initial value minus that consumed, 

(2-31) - AWf. (lbs) 

c rem ^ 

A time history comparison of the control fuel useage computed from 
measurements with that baaed on calculated total impulse and effective average 
specific impulse is output by the program. 
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2.2.4 Angular Accelerations 


Angular accelerations between control motor firings are used to compute 
the disturbing moments (equations 2-16 through 2-18). An optional input 
includes the slope and scale factors from a plotted time history or 
oscillograph record as shown in the sketch below. 



The angular acceleration is computed as, 

(2-32) 6 =■ Kp tan (deg/sec) 

where, 

IC(. is the paper speed in length units per second. 

Kp is the rate scale factor in degrees per second per length unit. 

X .0 is the slope of the rate trace in degrees. 
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3.0 PROGRAM DESCRIPTION 


3.1 General 


This computer program is programmed in FORTRAN IV for a GDC CYBER 175 
system. The coding is compatible with ANSI standards. It is arranged to 
operate with standard card input and line printer output. Optional CALCOMP 
plotting IS based on standard CALCOMP plotters and software. 

A main routine (UPSTAG) and four subroutines require approximately 19K 
words of computer memory. All output is stored in array variables to 
facilitate a well formatted print and CALCOMP plot output format. 

Program flow and user instructions are presented in the following 
paragraphs. Input and output of a sample problem is illustrated along with 
the detailed descriptions. 


3.2 Program Flow 

Program flow is straightforward in five basic parts: 

. input data 

. calculation of boost phase disturbances and impulse 
. calculation of coast phase variables 
. output data on line printer 
. optional CALCOMP plots 

A flow chart of the main routine (UPSTAG) is presented in Figure 4* A 
complete listing of the FORTRAN program and subroutines other than the 
standard CALCOMP library subroutines are presented in Appendix A. 

Descriptions of the four subroutines are presented in the following 
paragraphs. 


3.3 Subroutine Description 


Four subroutines are used to support the UPSTAG main program; C, CURVE, 
DASH, and TBLU. A brief description of each is presented below. 

£ 

This short subroutine computes the compressibility factor for dry 
diatomic nitrogen by a curve fit over the pressure range of 400 to 3000 psia 
and a temperature range of 50 to 100 degrees Fahrenheit. The call statement 
is. 
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CALL C(P, T, CR) 


where, 

P 13 the input nitrogen pressure (psia) 

T 13 the input nitrogen temperature (degrees Rankine) 

CR IS the output compressibility factor 

CURVE 

This subroutine is used to set up the sequence of CALCOMP plots for 
pitch, yaw, and roll disturbing moments and the effective thrust 
misalignments. Plots are set up for a 20 x 20 divisions per inch graph paper 
having a grid size of 10 i 7 l/2 inches (see the output data description in 
this section for a sample plot output). The data to be plotted is passed to 
this subroutine via common in arrays PVAR, YVAR, RVAR, and TVAH. Solid dashed 
and dash-dot line plots are accomplished by the subroutine DASH. The call 
statement is, 

CALL CURVE (NP, NQ, IROLL, NTITLE, DYP, DYR, DIET, DX) 

where, 

NP input number of time points in the arrays to be plotted 

NQ input controls 

NQ = 0 if aerodynamic moments are zero only the total 
disturbing moments are plotted 
NQ = 1 total disturbing moments in addition to aerodynamic 
and thrust misalignment moments are plotted 

IROLL input control integer 

IROLL = 0 roll disturbing moment not plotted 
IROLL = 1 roll disturbing moment is plotted 

NTITLE input title description of 80 characters with eight (8) 
'ten-letter' words 

DYP input ordinate scale factor for the pitch and yaw moments 

(ft-lbs per inch). 

DYR input ordinate scale factor for the roll moment (ft-lbs per 

inch) 

DYET input ordinate scale factor for thrust misalignment (degrees 

per inch) 

DT input abscissa scale factor for time (seconds per inch) 

Note that care must be taken in selecting scale factors so that plotted 
data falls on grid. Limiting of the plotted data is automatically invoked in 
CURVE through the call statements to DASH. 

DASH 


This subroutine plots a curve on a CALCOMP plotter for a set of ordinates 
and abscissas. The style and type of line drawn is selected by the user. 

Note that the CALCOMP plot is specified in inches; plotting on metric paper 
requires appropriate scaling change before entering this subroutine. 
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The call statement is. 


CALL DASH (X, Y, HP, Zl, Z2, SPACE, XSCALE, YSCALE, LSYMB, XLIM, YLIM) 


where, 

X 

Y 

NP 

Zl 

Z2 

SPACE 


XSCALE 

YSCALE 

LSYMB 


XLIM 

YLIM 


- input array of abscissa values 

- input array of ordinate values 

- number of points in X and Y to be plotted 

- for dashed-dot lines this is length of long line measured 
in inches (see sketch below) 

- for dashed-dot lines this is length of short line measured 
in inches (see sketch) 

- for dashed style lines this is the length of the space between 
lines measured in inches. 

SPACE = 0 gives a solid line plot 

SPACE = negative gives special CALCOMP symbols at each point 

- abscissa plot scale factor (units per inch) 

- ordinate plot scale factor (units per inch) 

- special CALCOMP symbol code number used if SPACE is negative 
(see code below) 

- (+) LSYMB gives straight solid lines between symbol points 

- (-) LSYMB gives only symbols at each point without lines 

- plot limiting of the abscissa (inches) points out of range, 
range will appear at this limit 

- plot range of ordinate (inches) 


For ease in use, the following styles are typically possible. 




LINE 


TYPE 

Zl 

Z2 

SPACE 

LSYMB 





Solid 

— 

-- 

0. 

0. 





Dashed 

0.25 

0.25 

0.10 

0. 

- - 

- - 

- - 

- - 

Dashed 

0.07 

0.07 

0.07 

0. 





Dashed 

0.5 

0.03 

0.07 

0. 





Dot 





/> 

A 


,A 

Symbols 



-0 1 

+2 

TU 

•“TTT 






A 

A 

A 

A 

Symbols 
(no line) 

— 

— 

-0.1 

-2 


TBLU 

This is a single table lookup subroutine. It is based on linear 
interpolation between points for a single array having alternating values of 
abscissas and ordinates. The abscissas must be in ascending order. 


The call to this subroutine is: 
CALL TBLU (NT, Y, X, T, M) 



NT - number of values in table 'T' including abscissas and 

ordinates. 

Y - 13 the ordinate to be found. 

X - IS the given abscissa. 

T - is the table of alternating abscissas and ordinates. 

M - is the table locator to begin the table search. M must be 

an integer value from 1 to NT. This index is changed by the 
subroutine to the location found in the lookup. 


3.4 Input Data Description 


Input data descriptions are presented in the following subparagraphs. A 
sample problem input data listing is presented in Figure 5 for reference. The 
input data can be separated into the following groups: 

1 ) option control card and title 

2 ) constants 

3 ) trajectoiy and vehicle characteristics tables 

4 ) angular accelerations (not shown in Figure 5) 

5 ) rate trace scale factors and slopes 

6) additional impulse table and incremental rates 

7 ) control fuel 

a) coast phase only constants (not shown in Figure 5) 

9 ) CALCOMP plot variables 


3 . 4.1 Options and Title 


The first card of input contains six integer control constants for the 
run. These are input in fields of five (5) columns. The number must be right 
justified. The input format is (lOI5)» 


FORTRAN 

NAME COLUMN 

lOPTl 5 


I0PT2 10 


DESCRIPTION 


option for input of pitch, yaw, and roll 
disturbing accelerations 
lOPTl = 0 input angular acceleration 
tables 

lOPTl = 1 input angular rate trace slopes 
and scale factors 

option for input of control fuel useage 
I0PT2 = 0 read in table of fuel remaining 
time history 

lOPTl = 1 read constants, pressure and 
temperature time histories and 
compute control fuel useage 
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FORTRAN 


NMB 

COLUMN 

DESCRIPTION 

I0PT3 

15 

coast input impulse option 
I0PT3 = 0 no incremental rates input 
I0PT3 = 1 input pitch yaw and roll incre- 
mental rates versus time during 
coast to compute impulse 

I0PT4 

20 

CALCOMP plot option 

I0PT4 = 0 no plot 

I0PT4 = 1 CALCOMP plot output 

I0PT5 

25 

boost and coast phase calculations 
I0PT5 = 0 boost phase only 
I0PT5 = •'■1 boost and coast phases 
I0PT5 = -1 coast phase only 

I0PT6 

30 

boost phase impulse option 
I0PT6 = 0 boost phase impulse computed 
from disturbing moments and 
additional impulse 

I0PT6 = 1 boost phase impulse computed 


from incremental rates 

The second and third card of input contains two lines of an arbitrary 
title which is output at the top of the printed page. Each card contains 80 
columns of hollerith data. It is read with a format of (SAIO). 


3 . 4.2 Constants 


Constants are input with fields of 10 columns. The fourth card contains 
five ( 5 ) time parameters as described below. 

FORTRAN 


NAME 

COLUMN 

DESCRIPTION 

TP 

sec 

time of stage ignition (flight time) 

TBO 

sec 

end of boost phase for stage being 
studied (flight time) 

TSTEP 

sec 

output step size desired during the boost 
phase (this is also the integration step 
size for impulse expended) 

T COAST 

sec 

time of coast phase termination (flight 
time) 

TSTEPC 

sec 

output step size for the coast phase 


12 


The next two cards of boost phase constants are input only if I0PT5 is 
2 sero or positive. The constants are read with format (6E10.3). 

FORTRAN 


NAME 

SYMBOL 

UNITS 

DESCRIPTION 

XT 

If 

inches 

body station of booster nozzle 
throat 

XC 

^c 

inches 

body station of pitch and yaw 
control motors 

ZC 

Zc 

inches 

location of pitch and yaw control 
motors away from centerline 

RC 

12 Rc 

inches 

roll control motor moment arm 

WH202I 

'^c 

lbs 

weight of control fuel onboard at 
stage ignition 

CNAS 

On s 

ft^/deg 

aerodynamic normal force 
coefficient slope at zero angle 
of attack times reference area 

XCP 

icp 

inches 

body station of aerodynamic 
center at zero angle of attack 

ETA 


degrees 

pitch and yaw control motor 
forward cant angle 

CN3 

a^S 

ft^/deg5 

cubic coefficient in aerodynamic 
normal force times reference area 

DCP 

9icp/3® 

inch/deg 

incremental change in aerodynamic 
center per degree angle of attack 


3 . 4.3 Boost Phase Tables (I0PT5 = 0, +l) 


When the boost phase disturbances are analyzed (I0PT5 = 0 or +l) the 
following eight (8) tables are input; dynamic pressure, angle of attack, angle 
of sideslip, booster thrust, weight of booster propellant remaining versus 
time, and roll moment of inertia, pitch or yaw moment of inertia, and center 
of mass versus percent propellant consumed. The first three tables of 
aerodynamic parameters are input with format (I5/(7E10.3) ) . The first card of 
each table is an integer number of values to be read into the table. This 
number is input in columns 1 through 5 and must be right justified. The table 
is entered with alternating values of abscissa and ordinate with seven (7) 
numbers per card. These three tables are dimensioned for 300 numbers. See 
Figure 5 for the sample problem input. The descriptions follow. 
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FORTRAN 

NAME 

NO. OF 
VALUES 

ABSCISSA 

UNITS 

ORDINATE 

UNITS 

DESCRIPTION 

NQ 

1 

— 

— 

— 

— 

number of values in 
table Q 

Q 

NQ 


sec 

Q 

Ibs/ft2 

dynamic pressure 
versus flight time 

NALPH 

1 

— 

— 

— 

- 

number of values in 
table ALPH 

ALPH 

NALPH 


sec 

a 

deg 

angle of attack ver: 
flight time 

NBETA 

1 

— 

— 

— 

- 

number of values in 
table BETA 

SETA 

NBETA 

tf 

sec 


deg 

angle of sideslip 
versus flight time 


The remaining input tables are dimensioned for 200 numbers each. These 
are read with format (I5/(6F10.3) ) • The first card of each table contains an 
integer number (right justified in columns 1 through 5) indicating the number 
of numbers included in the table. The table is read in with alternating 
abscissas and ordinates in fields of ten (lO) with siz (6) numbers per card. 
The descriptions of these tables follow . 


FORTRAN 

NAME 

NO. OF 
VALUES 

ABSCISSA 

UNITS 

ORDINATE 

UNITS 

DESCRIPTION 

NTT 

1 

— 

— 

— 

— 

number of values in 
table TT 

TT 

NTT 

ts 

sec 

T 

lbs 

booster thrust versus 
time after stage 
ignition 

NN¥P 

1 

— 

— 

— 

- 

number of values in 
table WP 

■rfP 

NNWP 


sec 

W 

Prem 

lbs 

booster propellant 
weight remaining 
versus time after 
booster ignition 

N¥X 

1 

— 

— 

— 

- 

number of values in 
table AIXX 

AIXX 

NWX 

^^cons 



slug-ft^ 

roll moment of inertia 
versus percent of 
booster propellant 
consumed 

NWY 

1 

— 

— 

— 

- 

number of values in 
table AIYY 
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FORTRAN NO. OF 


NAME 

VALUES 

ABSCISSA UNITS 

ORDINATE 

UNITS 

DESCRIPTION 

AIYY 

NWY 

^^cons 


slug-ft^ 

pitch or yaw moment of 
inertia versus percent 
of booster propellant 
consumed 

NXCG 

1 

— 

— 

- 

number of values in 
table XCG 

XCG 

NXCG 

^^cons 

^Gg 

inches 

center of mass body 
station versus percent 
of booster propellant 
consumed 


3.4*4 Boost Phase Angular Accelerations (lOPTl = O) 


Angular acceleration tables are input only when the boost phase is 
analyzed (I0PT5 - 0 or +l) and the acceleration input option lOPTl = 0. These 
are read in with format (I5/(6E10.3) )• 

FORTRAN NO. OP 


NAME 

VALUES 

ABSCISSA 

UNITS 

ORDINATE 

UNITS 

DESCRIPTION 

NTHE 

1 

— 

— 

— 

- 

number of values in 
table THEDD 

THEDD 

NTHE 

^f 

sec 

«• 

d 

deg/sec^ 

pitch angular 
acceleration versus 
flight time when 
control motors are off 

NPSI 

1 

— 

— 

— 

- 

number of values in 
table PSIDD 

PSIDD 

NPSI 

^f 

sec 


deg/sec^ 

yaw angular 
acceleration versus 
flight time when 
control motors are off 


NPHI 

1 


" 


number of values in 
table PHIDD 

PHIDD 

NPHI 


sec 

0 

deg/sec^ roll angular 

acceleration versus 
flight time when 
control motors are off 
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5.4.5 Boost Phase Rate Slopes (lOPTl = l) 


When lOPTl = 1 and I0PT5 = 0 or +1 the angular accelerations are computed 
fom the slopes of the rate traces and their scale factors (see paragraph 2.2.4 
and equation 2-32). A single card containing the six (6) rate trace scale 
factors 13 input before the slope table. These are input with format (6E10.3) 
and have the following definitions: 


FORTRAN 

NAIffl 

UNITS 

COLUMNS 

DESCRIPTION 

XKTHE 

deg/sec/inch 

1-10 

pitch rate trace scale factor 

XKPSI 

deg/ sec/inch 

11-20 

yaw rate trace scale factor 

XKPHI 

deg/sec/inch 

21-30 

roll rate trace scale factor 

XKTP 

inches/sec 

31-40 

pitch rate trace paper speed 
(inverse of time scale 
factor) 

XKTY 

inches/sec 

41-50 

yaw rate trace paper speed 
(inverse of time scale 
factor) 

XKTR 

inches/sec 

51-60 

roll rate trace paper speed 
(inverse of time scale 
factor) 


Three tables of the pitch, yaw and roll rate trace slopes are entered next 
with format (I5/(6E10.3) ) • 


FORTRAN 

NAME 

NO. OF 
VALUES 

ABSCISSA 

UNITS 

ORDINATE 

UNITS 

DESCRIPTION 

NTHE 

1 

— 

— 

— 

- 

number of values in 
THEDD 

THEDD 

NTHE 

tf 

sec 


deg 

pitch rate trace slope 
versus flight time when 
control motors are off 

NPSI 

1 

— 

— 

— 

- 

number of values in 
table PSIDD 

PSIDD 

NPSI 

tf 

sec 


deg 

yaw rate trace slope 
versus flight time 
when control motors 
are off 

NPHI 

1 

— 

— 

— 

- 

number of values in 
table PHIDD 

PHIDD 

NPHI 

tf 

sec 


deg 

roll rate trace slope 
versus flight time 
when control motors 
are off 
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3.4*6 Incremental Impulse 


Total impulse calculations include allowance for an incremental impulse 
which is precalculated. This allows for impulse not covered by moment balance 
during boost in any of three axes of control. It is up to the user to define 
entries into this table. This table is input with format (I5/6E10.3). The 
impulse is cimulative or a continuous sum. 

FORTRAN HO. OF 


NAME 

VALUES 

ABSCISSA 

UNITS 

ORDINATE 

UNITS 

DESCRIPTION 

NRI 

1 

— 

— 

— 

- 

number of values in 
table RI 

RI 

NRI 

tf 

sec 

It 

lb-sec 

sum of incremental 
impulse versus flight 
time 


When (I0PT3 = l) incremental rates are also used to compute impulse, such 
as, during limit cycle motion. Additional tables are required. These are 
input as an incremental absolute value of rate changes between time points 
(not cumulative) * These incremental rates are summed in the program to 
estimate total impulse. 

FORTRAN NO. OF 


NAME 

VALUES 

ABSCISSA 

UNITS 

ORDINATE 

UNITS 

DESCRIPTION 

NDPR 

1 

— 

— 

— 

- 

number of values in 
table DPR 

DPR 

NDPR 

tf 

sec 


deg/sec 

incremental pitch rate 
between time points 
versus flight time 
(absolute values) 

NDYR 

1 

— 

— 

— 

- 

number of values in 
table DYR 

DYR 

NDYR 

^f 

sec 


deg/sec 

incremental yaw rate 


between time points 
versus flight time 
(absolute values) 


NDRR 

1 




■ 

number of values in 
table DRR 

DRR 

NDRR 

tf 

sec 

A0 

deg/sec 

incremental roll rate 
between time points 
versus flight time 
(absolute values) 
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3.4.7 Control Fuel 


When I0PT2 = 0 the control fuel remaining versus flight time is input with 
format (I5/(6E10.3)) • 

FORTRAN NO. OF 


NAME 

VALUES 

ABSCISSA 

UNITS 

ORDINATE 

UNITS 

DESCRIPTION 

NW 

1 

— 

— 

— 

- 

number of values in 
table W 

W 

NW 

^f 

sec 

Wc 

lbs 

control fuel weight 
remaining versus 
flight time 


When I0PT2 = 1 the control fuel remaining is computed from the unregulated 
pressure and temperature of the pressurizing gas. The alternate input 
includes a card containing five (5) constants followed by two tables. These 
are, 

FORTRAN 


NAME 

COLUMN 

SYMBOL 

UNITS 

DESCRIPTION 

DSNS 

1-10 

Pl 

Ibs/in^ 

density of liquid control 
fuel 

COMFAC 

11-20 

^reg 

— 

compressibility factor of 
gas on regulated side 

PRSREG 

21-30 

^reg 

psia 

regulated pressure of system 

DVOL 

31-40 

^^reg 

in^ 

initial volume of gas on 
regulated side of system 

V0LN2 

41-50 


in^ 

volume of unregulated side 
of pressurization system 


The following tables are then entered with format (I5/6E10.3)) • 
FORTRAN NO. OP 

NAME VALUiS ABSCISSA UNITS ORDINATE UNITS DESCRIPTION 


NNP 1 

P NNP tf sec Pg 

NNT 1 -- 

T NNT tf sec Tg 


number of values in 
table P 

psia unregulated gas 

pressure versus flight 
time 

number of values in 
table T 

°F unregulated gas 

temperature versus 
flight time 
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3.4.8 Coast Phase Constants (I0PT5 = -l) 


When no "boost phase is analyzed (I0PT5 = -l) nine (9) constants are input 
with format (6E10.3). These are. 


FORTRAN 

NAME 

COLUMN 

SYMBOL 

UNITS 

DESCRIPTION 

RC 

1-10 

12Rc 

inches 

roll control motor moment 
arm 

WH202I 

11-20 

o 

o 

lbs 

value of control fuel 
onboard at stage ignition 

YMA 

21-30 

12 -Ic 

inches 

pitch or yaw control motor 
moment aim 

RIXX 

31-40 

Ix 

slug-ft^ 

roll moment of inertia 

YIYY 

41-50 


slug-ft^ 

pitch or yaw moment of 
inertia 

WH202M 

51-60 

Vc 

lbs 

control fuel onboard at 
stage burnout 

(next card) 





WH202C 

1-10 


lbs 

calculated control fuel at 
stage burnout 

BOIMP 

11-20 

It 

lb-sec 

total impulse expended 
during boost phase 

RIC 

21-30 

It 

lb-sec 

incremental impulse 
expended during the boost 
phase 


3.4.9 GALCOMP Plot Variables (I0PT4 = l) 


This group of input data is required if the optional CALCOMP plots are 
desired (I0PT4 = l). It includes scale factors, plot options and a title. 
The first card in this group is the title card containing 80 columns of 
alphanumeric infonnation to be included at the top of each plot. It is read 
into array NTITLE with format (8A10). The next card contains the variable 
IROLL in column 5. This is an integer control constant; 

IROLL = 0 no roll moment plotted 

IROLL = 1 roll moment time history is plotted 
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The third card contains four scale factor constants. These are read with 
format (4E10.3) and are, 

FORTRAN 


NAME 

COLUMNS 

UNITS 

DESCRIPTION 

SFT 

1-10 

sec/in 

abscissa time scale factor 
of plots (note plot is 
limited to 7 inches along 
abscissa axis) 

SFPYM 

11-20 

ft-lbs/in 

pitch and yaw moment scale 
factor (plot is limited to 
3 inches) 

SFET 

21-30 

deg/in 

thrust misalignment scale 
factor (plot is limited to 
+ 3 inches) 

SFRM 

31-40 

ft-lbs/in 

roll moment scale factor 
(plot 13 limited to _+ 5 
inches) 


Note that there is no plotted information if the boost phase is deleted 
(I0PT5 = -1). 


5*5 Output Data Description 


Output includes printed data and optional CALCOMP plots (if I0PT4 = l) • A 
detailed description of the output is presented in the following paragraphs 
with a sample problem output of Figures 6 and 7 for reference. The printed 
output occurs in five basic pages, 

. pitch boost phase variables 
. yaw boost phase variables 
. roll boost phase variables 
• boost phase system summary data 
. coast phase data 


3* 5.1 Boost Phase - Pitch (I0PT5 = 0, +l) 


This page of output includes the time histories of pitch variables during 
the boost phase. Both stage time (measured from ignition) and flight time 
(measured from liftoff) are included. Refer to Figure 6 for the output 
layout. The definitions of the printed names follow. 
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LABEL 

SYMBOL 

UNITS 

DESCRIPTION 

PLIGHT TIME 

tf 

sec 

flight time 

STAGE TIME 

^s 

sec 

stage time (measured from 
booster ignition) 

ANGULAR ACCEL 

e 

deg/sec^ 

pitch angular acceleration 

MOMENT OF INER 

h 

slug-ft^ 

pitch moment of inertia 

AERO MOMENT 

^aero 

ft-lbs 

aerodynamic pitching moment 
(Equation 2-ll) 

MISALIGN MOMENT 

Mr 

ft -lbs 

pitching moment due to 
booster (Equation 2-19) 

TOTAL MOMENT 

Mj) 

ft-lbs 

total pitch disturbing 
moment (Equation 2-16) 

THRUST MISALN 


deg 

pitch component of effective 
thrust misalignment 
(Equation 2-2l) 

IMPULSE 

It 

lb-sec 

total impulse due to pitch 
control motors during boost 
phase 


3»5«2 Boost 

Phase - Yaw (IOPT5 = 0, +l) 

This page format contains the following 

variables. 

LABEL 

SIDIBOL 

UNITS 

DESCRIPTION 

FLIGHT TIME 

tf 

sec 

flight time 

STAGE TIME 

“^s 

sec 

stage time measured from 
booster ignition 

ANGULAR ACCEL 


deg/sec^ 

yaw angular acceleration 

MOMENT OF INER 

h 

slug-ft^ 

pitch or yaw moment of 
inertia 

AERO MOMENT 

Maero 

ft-lbs 

yaw component of aerodynamic 
moment (Equation 2-12) 

MISALIGN MOMENT 

Nr 

ft-lbs 

yaw moment due to booster 
(Equation 2-20) 
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LABEL 

SYMBOL 

UNITS 

DESCRIPTION 

TOTAL MOMENT 

Nd 

ft-lbs 

total yaw disturbing moment 
(Equation 2-17) 

THRUST MISALN 


deg 

yaw component of effective 
thrust misalignment 
(Equation 2-22) 

LMPULSE 

It 

lb-sec 

yaw contribution to total 
impulse during boost phase 


3«5»5 Boost 

Phase - Roll (I0PT5 =0, +l) 

This page 

format includes the following variables, 

LABEL 

SYMBOL 

UNITS 

DESCRIPTION 

PLIGHT TIME 

^=f 

sec 

flight time 

STAGE TIME 

^3 

sec 

stage time measured from 
booster ignition 

ANGULAR ACCEL 


deg/sec^ 

roll angular acceleration 

MOMENT OF INER 

Ix 

slug-ft^ 

roll moment of inertia 

TOTAL MOMENT 

^D 

ft-lbs 

total roll disturbing moment 
(Equation 2-18) 

IMPULSE 

It 

lb-sec 

yaw contribution to total 
impulse during boost phase 

C.G. POINT 

^cg 

inches 

body station of center of 
mass 

CALCULATED 

CONSUM 

Wc 

lbs 

calculated control fuel 
consumed based on total 
pitch yaw and roll inpulse 
and average specific impulse 
(Equations 2-26) 

CALCULATED 

REMAIN 

w 

"c rem 

lbs 

calculated control fuel 
remaining based on total 
impulse computed and average 
specific impulse 
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3«5*4 Boost Phase System Variables (I0PT5 * 0, +l) 


This page 

format includes 

the following 

summary data. 

LABEL 

SIMBOL 

UNITS 

DESCRIPTION 

FLIGHT TIME 

tf 

sec 

flight time 

STAGE TIME 

ts 

sec 

stage time measured from 
booster ignition 

THRUST 

T 

lbs 

booster thrust 

TOTAL MISALIGN 

fr 

deg 

total effective booster 
thrust misalignment 

LAMDA 


deg 

roll orientation of thrust 
misalignment (Equation 2-24) 

SPECIFIC 

IMPULSE 

^sp 

sec 

effective specific impulse 
during the boost phase 
(Equation 2-28) 

TOTAL IMPULSE 

It 

lb-sec 

total pitch yaw and roll 
impulse expended (Equation 
2-26) 

FUEL CONSUM 

^c 

lbs 

control fuel cons\xmed 
measured or computed by 
Equation 2-29 

FUEL REMAIN 

^rem 

lbs 

control fuel remaining based 
on measured or Equation 2-29 


3*5*5 Coast Phase (I0PT5 “ + l) 


This page of output includes the coast phase performance and total flight 


results on total impulse. 



LABEL 

SYMBOL 

UNITS 

DESCRIPTION 

FLIGHT TIME 

^f 

sec 

flight time 

STAGE TIME 

^3 

sec 

stage time measured from 
booster ignition 

COAST IMPULSE 

It 

lb-sec 

total impulse expended 
during coast 
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LABEL 

SYMBOL 

UNITS 

DESCRIPTION 

TOTAL IMPULSE 

It 

lb-sec 

total impulse including 
boost and coast phases 

COAST ISP 

^sp 

Ibf- sec/lbjjj 

effective specific impulse 
during the coast phase 

TOTAL ISP 

^sp 

Ibf-sec/lbjj 

effective specific impulse 
including boost and coast 
phases 

COAST FUEL 

Wc 

lbs 

control fuel consumed during 
the coast phase only 

TOTAL FUEL 

^c 

lbs 

control fuel consumed 
including the boost and 
coast phases 

FUEL REM MEAS 

'"^rem 

lbs 

control fuel remaining based 
on measured data or Equation 
(2-29) 

FUEL REM CALC 

^rem 

lbs 

calculated control fuel 
remaining based on effective 
average specific impulse and 
impulse expended 


3.5.6 CALCOMP Plots (I0PT4 = l) 


The optional CALCOMP plots include the pitch, yaw and roll disturbing 
moments and the effective thrust misalignment versus flight time. Plots from 
the sample problem are presented in Figure 7. 
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Figure 1 
Sign Convention 
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Figure 3 

Typical Rate Time Histories 
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Figure 4 

Flow Chart of UPSTAG 
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Figure 4 (continued) 
Flow Chart of UPSTAG 
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Figure 4 (continued) 
Flow Chart of UPSTAG 
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Figure 4 (continued) 
Flow Chart of UPSTAG 
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Figure 4 (continued) 
Flow Chart of UPSTAG 
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Figure 4 (continued) 
Flow Chart of UPSTAG 
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Figure 5 

Sample Problem Inout Data 


1 1 0 

1 1 

1 0 



1 


SAnPLE PROBLEH FOR UPSTAG PROGRAM 




SCOUT UEHI 

CLE S-192C 

SECOND STAGE DISTURBING HONENTS 


36.16 

126.44 

0.25 

160.4 

1.0 



448.51 

467.76 

14.93 

16.89 

119.3 

0.472 


109.0 

5.0 

0.00358 

6.3 




322 


Q 

DYNAMIC 

PRESSURE 



84.000 

138.680 

84.500 

131.201 

85.000 

124.079 

85.500 

118.106 

86.000 

112.148 

86.500 

106.786 

87.000 

101.354 

37.500 

95.946 

88.000 

90.751 

88.500 

86.363 

89.000 

84.740 

39.500 

82.326 

90.000 

79.853 

90.500 

77.682 

91.000 

75.215 

91.500 

72.610 

92.000 

69.868 

92.500 

67.020 

93.000 

64.742 

93.500 

63.296 

94.000 

61.789 

94.500 

60.204 

95.000 

58.548 

95.500 

56.802 

96.000 

54.926 

96.500 

52.932 

97.000 

51.202 

97.500 

50.740 

98.000 

50.252 

98.500 

49.688 

99.000 

48.966 

99.500 

48.175 

100.000 

47.271 

100.500 

46.257 

101.000 

45.146 

101.500 

43.869 

102.000 

42.476 

102.500 

41.000 

103.000 

39.436 

103.500 

37.732 

104.000 

35.870 

104.500 

34.135 

105.000 

33.763 

105.500 

33.372 

106.000 

32.956 

106.500 

32.511 

107.000 

31.957 

107.500 

31.278 

108.000 

30.516 

108.500 

29.649 

109.000 

28.698 

109.500 

27.650 

110.000 

26.515 

110.500 

25.259 

111.000 

23.909 

111.500 

22.088 

112.000 

20.100 

112.500 

17.952 

113.000 

15.659 

113.500 

13.171 

114.000 

10.499 

114.500 

10.009 

115.000 

9.824 

115.500 

9.585 

116.000 

9.311 

116.500 

9.003 

117.000 

8.664 

117.500 

8.285 

118.000 

7.894 

118.500 

7.476 

119.000 

7.023 

119.500 

6.524 

120.000 

6.320 

120.500 

6.120 

121.000 

5.900 

121.500 

5.601 

122.000 

5.336 

122.500 

4.992 

123.000 

4.632 

123.500 

4.223 

124.000 

3.807 

124.500 

3.368 

125.000 

3.158 

125.500 

2.943 

126.000 

2.739 

126.500 

2.532 

127.000 

2.325 

127.500 

2.121 

128.000 

1.918 

128.500 

1.716 

129.000 

1.515 

129.500 

1.324 

130.000 

1.235 

130.500 

1.146 

131.000 

1.058 

131.500 

.971 

132.000 

.883 

132.500 

.797 

133.000 

.710 

133.500 

.625 

134.000 

.540 

134.500 

.455 

135.000 

.414 

135.500 

.386 

136.000 

.357 

136.500 

.329 

137.000 

.301 

137.500 

.273 

138.000 

.245 

138.500 

.218 

139.000 

.190 



222 


ALPHA 




84.000 

.087 

84.500 

.136 

85.000 

.134 

85.500 

.123 

86.000 

.100 

86.500 

.231 

87.000 

.167 

87.500 

.258 

88.000 

.252 

88.500 

.245 

89.000 

.212 

89.500 

.291 

90.000 

. 197 

90.500 

.306 

91 .000 

.181 

91.500 

.217 

92.000 

.149 

92.500 

.035 

93.000 

.163 

93.500 

.097 

94.000 

.246 

94.500 

.154 

95.000 

.215 

95.500 

.292 

96.000 

.217 

96.500 

.312 

97.000 

.333 

97.500 

.323 

98.000 

.462 

98.500 

.494 

99.000 

.449 

99.500 

.539 

100.000 

.588 

100.500 

.504 

101.000 

.531 

101.500 

.535 

102.000 

.428 

102.500 

.397 

103.000 

.509 

103.500 

.373 

104.000 

.273 

104.500 

.264 

105.000 

.213 

105.500 

.164 

106.000 

.210 

106.500 

.274 

107.000 

.215 

107.500 

.339 

108.000 

.335 

108.500 

.231 

109.000 

.303 

109.500 

.205 

110.000 

.140 

110.500 

.183 

111.000 

.024 

111.500 

.002 

112.000 

.050 

112.500 

-.003 

113.000 

.224 

113.500 

.050 

114.000 

.176 

114.500 

.159 

115.000 

.091 

115.500 

.211 

116.000 

.097 

116.500 

.153 

117.000 

.246 

117.500 

.046 

118.000 

.210 

118.500 

.061 

119.000 

.101 

119.500 

.272 

120.000 

.174 

120.500 

.371 

121.000 

.297 

121.500 

.182 

122.000 

.447 

122.500 

.011 

123.000 

.092 

123.500 

.077 

124.000 

-.632 

124.500 

-.465 

125.000 

.032 

125.500 

“ • S78 

126.000 

-.295 

126.500 

-.065 

127.000 

-.691 

127.500 

-.167 

128.000 

-.280 

128.500 

-.445 

129.000 

-.004 

129.500 

-.235 

130.000 

-.480 

130.500 

-.753 

131.000 

-.518 

131.500 

-.077 

132.000 

-.215 

132.500 

-.411 

133.000 

-.636 

133.500 

-.854 

134.000 

-.397 

134.500 

-.074 

135.000 

-.246 

135.500 

-.418 

136.000 

-.597 

136.500 

-.774 

137.000 

-.929 

137.500 

-.291 

138.000 

-.136 

138.500 

-.316 

139.000 

-.437 
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Figure b iconti nuea; 
Samole Problem Input Data 


222 

84 . 000 

-.440 

84.500 -.417 

85.000 

-.359 

85.500 

- . 400 

86.000 

-.386 

86.500 

-.351 

87.000 

-.136 

87.500 

.052 

88.000 

-.023 

88.500 

.021 

89.000 

. 248 

89.500 

.147 

90.000 

.031 

90.500 

.185 

91 . 000 

.276 

91.500 

.013 

92.000 

-.161 

92.500 

-.279 

93.000 

-.273 

93.500 

-.216 

94.000 

-.083 

94.500 

-.107 

95.000 

-.284 

95.500 

-.341 

96.000 

-.315 

96.500 

-.262 

97.000 

-.140 

97.500 

-.161 

98.000 

-.318 

98.500 

-.335 

99.000 

-.306 

99.500 

-.138 

100.000 

.141 

100.500 

.075 

101.000 

-.048 

101.500 

-.067 

102.000 

-.011 

102.500 

.127 

103.000 

-.112 

103.500 

-.333 

104.000 

-.492 

104.500 

-.592 

105.000 

-.580 

105.500 

-.400 

106.000 

-.140 

106.500 

-.389 

107.000 

-.477 

107.500 

-.472 

108.000 

-.303 

108.500 

-.180 

109.000 

-.013 

109.500 

.137 

110.000 

. 337 

110.500 

.238 

111.000 

-.107 

111.500 

- . 557 

112.000 

-.831 

112.500 

-.587 

113.000 

-.398 

113.500 

-.224 

114.000 

.044 

114.500 

.163 

115.000 

.219 

115.500 

-.429 

116.000 

-1.006 

116.500 

-.818 

117.000 

-.406 

117.500 

.022 

118.000 

-.013 

118.500 

-.440 

119.000 

- . 590 

119.500 

-.055 

120.000 

.505 

120.500 

.304 

121.000 

-.083 

121.500 

.109 

122.000 

. 886 

122.500 

. 151 

123.000 

-.090 

123.500 

.094 

124.000 

-.643 

124.500 

-.292 

125.000 

-.158 

125.500 

- . 521 

126.000 

.069 

126.500 

.092 

127.000 

-.322 

127.500 

.292 

128.000 

.135 

128.500 

-.310 

129.000 

-.088 

129.500 

.284 

130.000 

.459 

130.500 

.065 

131.000 

-.348 

131.500 

.014 

132.000 

.450 

132.500 

.339 

133.000 

.062 

133.500 

-.339 

134.000 

-.013 

134.500 

.400 

135.000 

.227 

135.500 

-.094 

136.000 

-.433 

136.500 

-.081 

137.000 

.437 

137.500 

.100 

138.000 

-.206 

138.500 

-.133 

139.000 

.404 



68 

BOOSTER 

THRUST UERSUS TIRE ( 

CASTOR II) 


0.00 

0.0 

.21 

42713.6 

.46 

40156.7 


1 . 86 

44017.7 

3.36 

46204.6 

6.36 

51153.7 


8 . 36 

54410.5 

10.36 

57460.7 

12.36 

60193.0 


14.36 

62885.3 

16.36 

65194.0 

18.36 

67299.8 


20.36 

68908.8 

22.36 

70342.4 

23.86 

71204.5 


25.36 

71801.8 

27.86 

71612.4 

29.36 

71214.3 


30 . 86 

70392.3 

31.86 

69206.1 

34.96 

69507.8 


35.56 

67122.7 

35.86 

66835.9 

35.96 

65433.8 


36.46 

27492.7 

36.66 

15760.3 

36.86 

10049.3 


37.26 

5339.6 

37.86 

1697.8 

38.36 

440.0 


39.41 

21.0 

40.00 

0.0 

100.00 

0.0 


66 

BOOSTER 

UEIGHT OF 

PROPELLANT 

RERAINING TIRE HISTORY 

0.00 

3274.34 

.21 

8253.70 

.46 

8213.68 


1.86 

7984.04 

3.36 

7733.81 

6.36 

7197.29 


8.36 

6805.81 

10.36 

6396.28 

12,36 

5967.30 


14.36 

5511.28 

16.36 

5041.14 

18.36 

4557.46 


20.36 

4056.87 

22.36 

3547.13 

23.86 

3159.52 


25.86 

2643.79 

27.86 

2130.12 

29.36 

1746.67 


30.86 

1372.14 

31.86 

1126.46 

34.96 

372.10 


35.56 

227.09 

35.86 

157.85 

35.96 

133.41 


36.46 

47.56 

36.66 

31.58 

36.86 

22.04 


37.26 

10.67 

37.86 

2.87 

38.36 

.89 


39.41 

0.00 

40.00 

0.00 

100.00 

0.00 


10 

IXX ROLL 

INERTIA 

US. % PROPELLANT CONSURED 


0. 

427.80 

25. 

410.55 

50. 

367.84 


75. 

299.68 

100. 

206.10 




10 

IVY PITCH/YAU INERTIA US. X 

PROPELLANT CONSURED 


0 . 

49263.29 

25. 

45824.92 

50. 

41600.69 


75. 

36058.32 

100. 

28031.35 





1® XCG us. % PROPELLANT CONSUPIED 

0. 382.03 25. 272.09 50. 258.44 

75. 238.53 100. 206.76 



Figure 5 (concluded) 
Sample Problem Inout Data 


2.82S 

3.841 

3.318 

3.96 

3.96 

3.96 


1 OS 


PITCH SLOPES 




8S.36 

5.90 

86.76 

5.00 

87,26 

9.40 


S7.8fi 

7.60 

88.86 

6.30 

89.86 

6.90 


90.68 

7.50 

91.16 

8.40 

91.76 

8.90 


93.76 

8.90 

93.56 

9.50 

94.46 

9.50 


95.36 

10.50 

95.96 

11.10 

96.86 

11.50 


97,56 

11.50 

98,36 

11.50 

99.06 

12,50 


99.86 

13.30 

100.46 

13.50 

101.36 

13.00 

— 

101 .96 

13.70 

103.76 

13.00 

103.36 

13.80 


10'4.16 

13.90 

104.96 

13.00 

105.66 

13.10 


i06.46 

13.00 

107.16 

13.70 

107.86 

14.30 


108.66 

13.30 

109.36 

14.30 

110.16 

14.50 


110.96 

13.60 

111.86 

13,50 

113.66 

13.40 


113.46 

13.60 

114.38 

14.30 

115.16 

15.00 


116.46 

15,10 

117,06 

15.80 

117.56 

15.90 


118.46 

16.90 

119.36 

17.30 

130.06 

18.00 


130.96 

19.00 

131.46 

30.00 

131.76 

31.30 


133.36 

16.00 

133.66 

11.10 

133.06 

3.30 


133.66 

.80 

134.76 

0.00 

186.16 

0.00 


98 


YAU ■ 

SLOPES 




86.36 

-3.50 

87.16 

-4.90 

87.96 

-3.80 


88.96 

-4.90 

89.46 

-3.90 

90.16 

-1.90 


90.66 

-3.50 

91.46 

-3,00 

93.06 

-1.30 


93.16 

-1,00 

94.06 

-3.40 

94.66 

-3.60 


95.66 

-1.90 

95.36 

-3.10 

97.06 

-3.00 


97.46 

-3.10 

93.66 

-3.40 

99.66 

-3.10 


100.46 

-3.80 

101.16 

-3.50 

103.16 

-2.39 


103.96 

-3.30 

103.96 

-3.00 

105.16 

-1.80 


105.76 

-3.30 

106.56 

-3.00 

107.16 

-3.40 


108.16 

-1.00 

109.16 

-.30 

110.16 

-.80 


111.16 

0.00 

113.16 

.80 

113.16 

1.00 


114.16 

.30 

115.66 

.30 

116.66 

.30 


117.36 

.10 

119.36 

.50 

130.66 

* 69 


131.66 

-.30 

133.36 

3.50 

133.41 

9.00 


133.66 

-1.00 

133.06 

-2.70 

133.36 

.70 


133.86 

1.00 

134.46 

-.30 

135.16 

0.00 


186.16 

0.00 






53 


ROLL 

SLOPES 




86.96 

3.50 

87.96 

1.30 

88.76 

* 19 


89.96 

1.10 

91.66 

1.00 

93.56 

1.00 


95.16 

1.50 

96.76 

1.10 

98.36 

1.00 


99.06 

1.50 

99.76 

1,30 

101.16 

1.00 


103.36 

1.30 

104.86 

.90 

105.46 

0.00 


107.06 

-.60 

107.76 

-1,00 

108.56 

-1.40 


109.16 

-1.40 

110.86 

-.60 

111.76 

-2.00 


113.76 

-.60 

113.46 

-.40 

114.16 

0.00 


118.16 

0.00 

131.66 

0.00 

133.36 

3.00 


133.16 

0.00 

186.16 

0.00 




33 

INCREnEHTAL IMPULSE 

US. FLIGHT TINE 



86 . 

0.0 

98. 

0.0 

103. 

3.8 


106. 

3.9 

110. 

7.3 

114. 

63.0 


118. 

309.3 

133. 

419.0 

134. 

530.0 


136. 

931.0 

160. 

931.0 




33 

PITCH DELTA RATES DURING COAST 



0. 

0.00 

136. 

0.00 

130. 

8.94 


134, 

5.44 

138. 

3.48 

143. 

3.56 


146. 

0.00 

150. 

3.87 

154. 

0.00 


156. 

5.09 

160. 

1.93 




33 

YAU DELTA RATES DURING COAST 




0. 

0.00 

136. 

0.00 

130. 

10,63 


134. 

5.11 

138. 

5.33 

143. 

7.37 


146. 

6.99 

150. 

5.33 

154. 

7.44 


156. 

5.34 

160. 

3.58 




4 

ROLL DELTA RA^ES DURING COAST 



0. 

0. 

160. 

0. 




33 

CONTROL 

FUEL RENAINING US. FLIGHT TINE 



80. 

119.4 

90. 

119.35 

95. 

114.7 


100. 

109.7 

105. 

104.3 

110. 

99.8 


115. 

95.8 

130. 

90.3 

125. 

84.6 


130. 

78.3 

135. 

75,3 

140. 

73.3 


145. 

70.1 

150. 

68.5 

IS5. 

57.6 


160. 

66.0 






SrtflPLE PROBLEfl - SCOUT 

S-193C SECOND STAGE 

DISTURBANCES 


1 

CALCOflP 

PLOT OPTION 

WITH ROLL 

. PLOT 



10. 

3000. 

0.1 

3. 




lEOR 









Figure 6 

Sample Problem Printed Output Data 


UPPER STAGE MOMENT ROUTINE PAGE - 1 A 

PITCH UARIABLES TIME HISTORIES 
SAMPLE PROBLEM FOR UPSTAG PROGRAM 

SCOUT UEHICLE S-193C SECOND STAGE DISTURBING MOMENTS 


FLIGHT 

STAGE 

ANGULAR 

MOMENT 

AERO 

MISALIGN 

TOTAL 

THRUST 


TIME 

TIME 

ACCEL 

OF INER. 

MOMENT 

MOMENT 

MOMENT 

MISALN 

IMPULSE 

SEC 

SEC 

DEG^S/'S 

SLUG-FTa 

FT-LB 

FT-LB 

FT-LB 

DEG 

LB-SEC 


8G.16 

0.00 

1.16 

49263.3 

105.23 

888.61 

993.84 

0.0000 

0.00 

87.00 

.84 

1.40 

49058.9 

113.94 

1082.82 

1196.76 

.1082 

59.07 

88.00 

1.84 

1.45 

48786.2 

153.04 

1084.68 

1237.73 

.1011 

136.89 

89.00 

2.84 

1.23 

48509.0 

119.37 

925.81 

1045.18 

.0831 

209.56 

90.00 

3.84 

1.37 

48222.1 

104.41 

1052.34 

1156.74 

.0908 

279.34 

91.00 

4.84 

1.59 

47924.9 

89.54 

1244.08 

1333.62 

.1032 

357.91 

98.00 

5.84 

1.75 

47627.6 

68.37 

1387.64 

1456.01 

.1108 

445.52 

93.00 

6.84 

1.79 

47316.9 

68.71 

1407.59 

1476.30 

.1083 

537.19 

94.00 

7.84 

1.87 

46991.5 

98.60 

1436.56 

1535.16 

.1066 

630.88 

95.00 

8.84 

2.01 

46659.0 

80.91 

1554.03 

1634.94 

.1114 

729.01 

96.00 

9.84 

2.20 

46318.6 

76.08 

1700.88 

1776.96 

.1179 

834.10 

97.00 

10.84 

2.28 

45970.5 

108.24 

1717.62 

1825.86 

.1155 

944.50 

98.00 

11.84 

2.28 

45565.7 

145.28 

1664.50 

1809.79 

.1085 

1055.19 

99.00 

12.84 

2.46 

45114.5 

136.39 

1802.36 

1938.76 

.1140 

1168.49 

100.00 

13.84 

2.43 

44648.9 

170.46 

1725.25 

1895.71 

.1058 

1283.52 

101.00 

14.84 

2.55 

44176.3 

145.90 

1819.42 

1965.33 

. 1085 

1398.47 

102.00 

15.84 

2.52 

43696.3 

109.93 

1814.84 

1924.76 

.1054 

1513.40 

103.00 

16.84 

2.57 

43209.6 

119.75 

1815.33 

1935.07 

. 1028 

1626.57 

104.00 

17.84 

2.56 

42715.8 

57.72 

1849.10 

1906.82 

.1021 

1738.34 

105.00 

18.84 

2.58 

42213.6 

41.85 

1861.59 

1903.44 

. 1005 

1848.33 

106.00 

19.84 

2.59 

41702.5 

40.36 

1847.79 

1883.15 

.0978 

1956.92 

107.00 

20.84 

2.69 

41057.8 

39.05 

1891.23 

1930.28 

.0978 

2065.11 

108.00 

21.84 

2.79 

40374.9 

57.26 

1911.69 

1968.95 

.0966 

2174.32 

109.00 

22.84 

2.73 

39687.5 

48.14 

1845.67 

1893.81 

.0912 

2281.29 

110.00 

23.84 

2.38 

38995.1 

20.14 

1942.96 

1963.10 

.0941 

2386.89 

111.00 

24.84 

2,50 

38304.2 

3.09 

1667.78 

1670.87 

.0794 

2485.27 

112.00 

25.84 

2,51 

37613.3 

5.16 

1643.97 

1649.13 

.0770 

2574.16 

113.00 

26.84 

2,68 

36925.1 

17.93 

1710.70 

1728.63 

.0793 

2663.62 

114.00 

27.84 

2.78 

36237.0 

9.38 

1749.22 

1758.60 

.0802 

2754.98 

115.00 

28.84 

2.96 

35325.0 

4.40 

1822.85 

1827.26 

.0825 

2847.49 

116.00 

29.84 

3.01 

34344.1 

4.17 

1800.45 

1804.62 

.0805 

2939.67 

117.00 

30.84 

3.15 

33375.2 

9.73 

1825.40 

1835.14 

.0808 

3030.57 

118.00 

31.84 

3.29 

32421.6 

7.42 

1854.23 

1861.65 

.0821 

3121.48 

119.00 

32.84 

3,44 

31477.1 

3.01 

1887.73 

1890.75 

.0821 

3212.36 

120.00 

33.84 

3.62 

30532.8 

4.52 

1925.33 

1929.85 

.0823 

3303.53 

121.00 

34.84 

3.87 

29588.6 

7.03 

1990.98 

1998.01 

.0836 

3395.89 

122.00 

35.84 

3.90 

28661.8 

8.88 

1941.81 

1950.70 

.0835 

3487.43 

123.00 

36.84 

.86 

28120.6 

1.63 

420.81 

422.44 

.1129 

3541.99 

124.00 

37.84 

.11 

28043.5 

-8.77 

61.59 

52.82 

.0963 

3552.91 

125.00 

38.84 

0.00 

28033.2 

.33 

-.38 

0.00 

-.0044 

3554.12 

126.00 

39.84 

0.00 

28031.4 

-3.05 

3.05 

0.00 

1.5225 

3554.12 

126.44 

40.28 

0.00 

28031.4 

-.90 

.90 

0.00 

0.0000 

3554.12 


38 



Figure 6 (continued) 

Sample Problem Printed Output Data 


UPPER STAGE nOPlENT ROUTINE PAGE -IB 

YAU UARIABLES TIME HISTORIES 
SAMPLE PROBLEM FOR UPSTAG PROGRAM 

SCOUT UEHICLE S-193C SECOND STAGE DISTURBING MOMENTS 


FLIGHT 

STAGE 

ANGULAR 

MOMENT 

AERO 

MISALIGN 

TOTAL 

THRUST 


TIME 

TIME 

ACCEL 

OF INER. 

MOMENT 

MOMENT 

MOMENT 

MISALN 

IMPULSE 

SEC 

SEC 

DEG/'S/S 

SLUG-FTE 

FT-LB 

FT-LB 

FT-LB 

DEG 

LB-SEC 


86.16 

0.00 

-.49 

49263.3 

277.92 

-700.20 

-422.28 

0.0000 

0.00 

87,00 

.84 

-.88 

49058.9 

92.79 

-846.47 

-753.68 

-.0845 

31.71 

88.00 

1.84 

-.57 

48786.2 

13.97 

-496.52 

-482.55 

-.0463 

71.23 

89.00 

2.84 

-.95 

48509.0 

-139.64 

-663.45 

-803.09 

-.0595 

112.15 

90.00 

3.84 

-.46 

48222.1 

-42.93 

-346.87 

-389.80 

-.0299 

149.95 

91.00 

4.84 

-.53 

47924.9 

-136.54 

-309.24 

-445.78 

-.0257 

176.32 

98.00 

5.84 

-.27 

47627.6 

73.88 

-299.16 

-225.28 

-.0239 

197.39 

93.00 

6.84 

-.20 

47316.9 

115.07 

-281.94 

-166.87 

-.0217 

209.65 

94.00 

7.84 

-.45 

46991.5 

33.27 

-404.89 

-371.63 

-.0300 

226.40 

95.00 

8.84 

-.46 

46659.0 

106.88 

-484.74 

-377.86 

-.0347 

249.61 

96.00 

9.84 

-.39 

46318.6 

110.44 

-427.54 

-317.10 

-.0296 

271.01 

97.00 

10.84 

-.57 

45970.5 

45.51 

-506.32 

-460.82 

-.0340 

294 . 85 

98.00 

11.84 

-.44 

45565.7 

100.00 

-449. 14 

-349. 14 

-.0293 

319.51 

99.00 

12.84 

-.45 

45114.5 

92.95 

-448.39 

-355.43 

-.0284 

340.80 

100.00 

13.84 

-.55 

44648.9 

-40.87 

-391.39 

—432 . 26 

-.0240 

364.43 

101.00 

14.84 

-.55 

44176.3 

13.19 

-436.98 

-423.79 

-.0261 

389.92 

102.00 

15.84 

-.44 

43696.3 

2.83 

-339.59 

-336.76 

-.0197 

412.39 

103.00 

16.84 

-.45 

43209.6 

26.35 

-365.29 

-338.94 

-.0207 

432.20 

10^.00 

17.84 

-.39 

42715.8 

104.02 

-395.90 

-291.88 

-.0219 

450.55 

105.00 

18.84 

-.36 

42213.6 

113.96 

-378.27 

-264.31 

-.0204 

466.61 

106.00 

19.84 

-.48 

41702.5 

26.90 

-375.80 

-348.89 

-.0199 

484.17 

107.00 

20.84 

-.50 

41057.8 

86.63 

-447.04 

-360.41 

-.0231 

504.27 

108.00 

21.84 

-.24 

40374.9 

51.79 

-221.17 

-169.38 

-.0112 

519.11 

109.00 

22.84 

-.08 . 

39687.5 

2.07 

-58.10 

-56.03 

-.0029 

525.35 

110.00 

23.84 

-.14 

38995.1 

-48.47 

-47.74 

-96.21 

-.0023 

529.52 

111.00 

24.84 

-.03 

38304.2 

13.78 

-30.58 

-16.80 

-.0015 

532.58 

112.00 

25.84 

.13 

37613.3 

85.77 

.84 

86.62 

.0000 

534.45 

113.00 

26.84 

.19 

36925.1 

31.86 

90.63 

122.49 

.0042 

539.98 

114.00 

27.84 

.08 

36237.0 

-2.35 

53.50 

51 . 16 

.0025 

544.53 

115.00 

28.84 

.05 

35325.0 

-10.59 

40.13 

29.53 

.0018 

546.62 

116.00 

29.84 

.04 

34344.1 

43.21 

-19.67 

23.54 

-.0009 

547.96 

117.00 

30.84 

.03 

33375.2 

16.07 

1.25 

17.32 

.0001 

548.98 

118.00 

31.84 

.04 

32421.6 

.46 

24.87 

25.33 

.0011 

550.03 

119.00 

32.84 

.08 

31477.1 

17.59 

28.58 

46.16 

.0012 

551.76 

120.00 

33.84 

.11 

30532.8 

-13.12 

70.58 

57.46 

.0030 

554.24 

121.00 

34.84 

.06 

29588.6 

1.97 

27.84 

29.81 

.0012 

556.29 

122.00 

35.84 

.36 

28661.8 

-17.61 

199.87 

182.26 

.0086 

561.20 

123.00 

36.84 

-.48 

28120.6 

1.60 

-237.35 

-235.76 

-.0637 

562.43 

124.00 

37.84 

.14 

28043.5 

8.92 

60.27 

69.19 

.0942 

566.26 

125.00 

38.84 

-.01 

28033.2 

1.90 

-6.29 

-4.39 

-.0720 

567.75 

126.00 

39.84 

0.00 

28031.4 

-.71 

.71 

0.00 

.3561 

567.85 

126.44 

40.28 

0.00 

28031.4 

-.87 

.87 

0.00 

0.0000 

567.85 
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Figure 6 (continued) 

Sample Problem Printed Output Data 


UPPER STAGE MOMENT ROUTINE PAGE -1C 

ROLL VARIABLES TIME HISTORIES 
SAMPLE PROBLEM FOR UPSTAG PROGRAM 

SCOUT VEHICLE S-192C SECOND STAGE DISTURBING MOMENTS 


FLIGHT 

STAGE 

ANGULAR 

MOMENT 

TOTAL 


C.G. 

CALCULATED FUEL 

TIME 

TIME 

ACCEL 

OF INER 

MOMENT 

IMPULSE 

POINT 

CONSUM 

REMAIN 

SEC 

SEC 

DEG/S/S 

SLUG-FT2 

FT-LB 

LB-SEC 

INCH 

LB 

LB 


86.16 

0.00 

.56 

427.8 

4.19 

0.00 

282.030 

0.000 

119.369 

87.00 

.84 

.55 

426.8 

4.07 

2.47 

281.439 

.669 

118.700 

88.00 

1.84 

.18 

425.4 

1.36 

1.93 

280.651 

1.507 

117.862 

89.00 

2.84 

.05 

424.0 

.36 

.61 

279.849 

2.313 

117.057 

90.00 

3.84 

.18 

422.6 

1.30 

.59 

279.020 

3.084 

116.285 

91.00 

4.84 

.17 

421.1 

1.22 

.90 

278.161 

3.839 

115.530 

92.00 

5.84 

.16 

419.6 

1.17 

.85 

277.301 

4.619 

114.750 

93.00 

6.84 

.16 

418.0 

1.17 

.83 

276.403 

5.364 

114.005 

94.00 

7.84 

.18 

416.4 

1.32 

.89 

275.462 

6.157 

113.212 

95.00 

8.84 

.23 

414.7 

1.68 

1.07 

274.501 

7.029 

112.340 

96.00 

9.84 

.21 

413.0 

1.49 

1.13 

273.517 

7.937 

111.432 

97.00 

10.84 

.17 

411.3 

1.25 

.97 

272.511 

8.899 

110.470 

98.00 

11.84 

.16 

407.9 

1.17 

.86 

271.252 

9.869 

109.500 

99.00 

12.84 

.23 

403.4 

1.64 

1.70 

269.794 

10.841 

108.528 

100.00 

13.84 

.19 

398.7 

1.30 

2.45 

268.290 

11.841 

107.528 

101.00 

14.84 

. 16 

393.9 

1.13 

2.96 

266.763 

12.852 

106.517 

102.00 

15.84 

.17 

389.0 

1.17 

3.62 

265.212 

13.843 

105.526 

103.00 

16.84 

.19 

384.1 

1.25 

3.94 

263.639 

14.799 

104.570 

104.00 

17.84 

.17 

379.1 

1.14 

4.20 

262.043 

15.735 

103.634 

105.00 

18.84 

.11 

374.0 

.72 

4.29 

260.421 

16.640 

102.729 

106.00 

19.84 

-.03 

368.9 

-.21 

4.08 

258.769 

17.543 

101.826 

107.00 

20.84 

-.09 

361.2 

-.58 

5.03 

256.490 

18.471 

100.899 

108.00 

21.84 

-.18 

352.8 

-1.10 

6.20 

254.037 

19.369 

100.000 

109.00 

22.84 

-.22 

344.3 

-1.35 

7.32 

251.567 

20.189 

99.180 

110.00 

23.84 

-.16 

335.8 

-.94 

8.11 

249.080 

20.983 

98.387 

111.00 

24.84 

-.13 

327.3 

-.75 

21.83 

246.598 

21.809 

97.560 

112.00 

25.84 

-.27 

318.3 

-1.48 

35.94 

244.116 

22.561 

96.808 

113.00 

26.84 

-.09 

310.3 

-.46 

49.77 

241.644 

23.342 

96.027 

114.00 

27.84 

-.01 

301.9 

-.08 

63.19 

239.172 

24.126 

95.243 

115.00 

28.84 

0.00 

291 . 1 

0.00 

99.60 

235.627 

25.066 

94.303 

116.00 

29.84 

0.00 

279.7 

0.00 

136.15 

231.745 

26.000 

93.370 

117.00 

30.84 

0.00 

268.4 

0.00 

172.73 

227.910 

26.922 

92.448 

118.00 

31.84 

0.00 

257.3 

0.00 

209.30 

224.136 

27.844 

91.525 

119.00 

32.84 

0.00 

246.3 

0.00 

261.73 

220.398 

28.884 

90.485 

120.00 

33.84 

0.00 

235.3 

0.00 

314.15 

216.661 

29.932 

89.437 

121.00 

34.84 

0.00 

224.3 

0.00 

366.58 

212.923 

30.986 

88.383 

122.00 

35.84 

. 16 

213.4 

.58 

419.21 

209.255 

32.056 

87.314 

123.00 

36.84 

.06 

207.1 

.23 

474.79 

207.113 

32.855 

86.515 

124.00 

37.84 

0.00 

206.2 

0.00 

530.08 

206.808 

33.357 

86.012 

125.00 

38.84 

0.00 

206.1 

0.00 

755.50 

206.767 

34.994 

84.375 

126.00 

39.84 

0.00 

206.1 

0.00 

981.00 

206.760 

36.612 

82 . 757 

126.44 

40.28 

0.00 

206.1 

0.00 

981.00 

206.760 

36.612 

82 . 757 
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Figure 6 (continued) 

Sample Problem Printed Output Data 


UPPER STAGE MOMENT ROUTINE PAGE -ID 

TIME HISTORIES OF SYSTEM 
SAMPLE PROBLEM FOR UPSTAG PROGRAM 

SCOUT UEHICLE S-ISEC SECOND STAGE DISTURBING MOMENTS 


FLIGHT 

STAGE 


TOTAL 


SPECIFIC 

TOTAL 

FUEL 

FUEL 

TIME 

TIME 

THRUST 

MISALGN 

LAMDA 

IMPULSE 

IMPULSE 

CONSUM 

REMAIN 

SEC 

SEC 

LB 

DEG 

DEG 

SEC 

LB-SEC 

LB 

LB 


86. 16 

0.00 

0 . 

0.0000 

0.00 

0.00 

0.00 

0.00 

119.37 

87.00 

.84 

41305 . 

.1373 

331.98 

xxxttxx 

93.34 

.00 

119.37 

88.00 

1.84 

43963 . 

.1113 

335.40 

xxxxxxx 

310.05 

.01 

119.36 

89.00 

3.84 

45446 , 

.1033 

334.37 

xxxxxxx 

333.33 

.01 

119.36 

90.00 

3.84 

46996 . 

.0956 

341.76 

xxxxxxx 

439.88 

.03 

119,35 

91.00 

4,84 

48646 . 

.1064 

346.04 

563.76 

535.13 

.95 

118.43 

93.00 

5.84 

50396 . 

.1133 

347.83 

343.58 

643.77 

1.88 

117.49 

93.00 

6.84 

51935 . 

.1104 

348,67 

366.15 

747.68 

3.81 

116.56 

94.00 

7.84 

53564 . 

.1107 

344.36 

339.51 

858.17 

3.74 

115.63 

95.00 

8.84 

55143 . 

.1167 

343.67 

309.83 

979.68 

4.67 

114.70 

96.00 

9.84 

56668 . 

.1316 

345.89 

195.13 

1106.33 

5.67 

113.70 

97.00 

10.84 

58118 . 

. 1304 

343.57 

185.98 

1340.33 

6.67 

113.70 

98.00 

11.84 

59486 . 

.1134 

344.90 

179.36 

1375.56 

7,67 

111.70 

99.00 

13.84 

60843 . 

.1174 

346.03 

174.39 

1510.99 

8.67 

110.70 

100.00 

13.84 

63187 . 

.1085 

347.33 

170.69 

1650.40 

9.67 

109.70 

101.00 

14.84 

63439 . 

.1116 

346.49 

166.34 

1791.35 

10.77 

108.60 

103.00 

15.84 

64594 . 

.1073 

349.40 

163.56 

1939.41 

11.87 

107.50 

103.00 

16.84 

65699 . 

.1048 

348.63 

159.05 

3063.70 

13.97 

106.40 

104.00 

17.84 

66753 . 

.1045 

347.91 

155.88 

3193.09 

14.07 

105.30 

105.00 

18.84 

67686 . 

.1036 

348.51 

153.89 

3319.33 

15.17 

104.30 

106.00 

19.84 

68490 . 

.0998 

348.50 

153.35 

3445.17 

16.05 

103.33 

107.00 

30.84 

69353 . 

.1005 

346.70 

153.07 

3574.41 

16.93 

103.44 

108.00 

31.84 

69970 . 

.0973 

353.40 

151.59 

3699.63 

17.81 

101.56 

109.00 

33.84 

70618 . 

.0913 

358.30 

150.57 

3813.96 

18.69 

100.68 

110.00 

33.84 

71193 . 

.0941 

358.59 

149.44 

3934.53 

19.57 

99.80 

111.00 

34.84 

71497 . 

.0795 

358.95 

149.33 

3039.67 

30.37 

99.00 

113.00 

35.84 

71796 . 

.0770 

.03 

148.54 

3144.55 

31.17 

98.30 

113.00 

36.84 

71709 . 

.0794 

3.03 

148.09 

3353.37 

31.97 

97.40 

114.00 

37.84 

71614 . 

.0803 

1.75 

147.69 

3363.70 

33.77 

96.60 

115.00 

38.84 

71353 . 

.0835 

1.36 

148.33 

3493.71 

33.57 

95.80 

116.00 

39.84 

70951 . 

.0805 

359.37 

146.89 

3633.78 

34.67 

94.70 

117.00 

30.84 

70403 . 

.0808 

.04 

145.61 

3753.38 

35.77 

93.60 

118.00 

31.84 

69330 . 

.0831 

.77 

144.43 

3880.81 

36.87 

93.50 

119.00 

33.84 

69301 . 

.0831 

.87 

143.94 

4035.85 

37.97 

91.40 

130,00 

33.84 

69399 . 

.0833 

3.10 

143.53 

4171.91 

39.07 

90.30 

131.00 

34.84 

69496 . 

.0836 

.80 

143.96 

4318.75 

30.31 

89.16 

133.00 

35.84 

66855 . 

.0839 

5.88 

143.53 

4467.84 

31.35 

88.03 

133.00 

36.84 

10630 . 

. 1396 

330.57 

140.95 

4579.31 

33.49 

86.88 

134,00 

37.84 

1819 . 

. 1348 

44.33 

138.35 

4649.35 

33.63 

85.74 

135.00 

38.84 

348 . 

.0731 

366.51 

140.38 

4877.37 

34.77 

84.60 

136.00 

39.84 

6 . 

1.5636 

13.17 

141.56 

5103.97 

36.05 

83.33 

136.44 

40.38 

0 . 

0.0000 

0.00 

139.38 

5103.97 

36.61 

83.76 


(Note: ******* results when the computed number overflows the available 
output field width) 
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Figure 6 (concluded) 

Sample Problem Printed Output Data 


UPPER STAGE flOHENT ROUTINE PAGE -IE 

COAST UARIABLES TIME HISTORIES 
SAMPLE PROBLEM FOR UPSTAG PROGRAM 

SCOUT UEHICLE S-19EC SECOND STAGE DISTURBING MOMENTS 


FLIGHT 

STAGE 

COAST 

TOTAL 

COAST 

TOTAL 

COAST 

TOTAL 

FUEL 

REMAINING 

TIME 

TIME 

IMPULSE 

IMPULSE 

ISP 

ISP 

FUEL 

FUEL 

MEAS 

CALC 

SEC 

SEC 

LB-SEC 

LB-SEC 

SEC 

SEC 

LB 

LB 

LB 

LB 


1E6.44 

40.38 

48.3 

5151.3 

0.00 

140.96 

0.00 

36.54 

83.76 

83.35 

137.00 

40.84 

109.9 

5313.8 

153.36 

139.90 

.73 

37.36 

83.04 

81.84 

139.00 

43.84 

339.6 

5433.5 

100.58 

136.43 

3.38 

39.83 

79.48 

80.00 

131.00 

44.84 

498.7 

5601.6 

99.01 

134.73 

5.04 

41 .58 

77.73 

78.58 

133.00 

46.84 

617.3 

5730.3 

103.93 

134.47 

6.00 

43.54 

76.76 

77.59 

135.00 

48.84 

735.3 

5838.3 

104.36 

133.98 

6.96 

43.50 

75.80 

76.68 

137.00 

50.84 

833.0 

5936.0 

103.44 

133.17 

7.96 

44.50 

74.80 

75.87 

139.00 

53.84 

933.7 

6035.7 

104.13 

133.65 

8.96 

45.50 

73.80 

74.95 

141.00 

54.84 

1054.4 

6157.3 

104.43 

133.03 

10.10 

46.64 

^ 73.66 

73.93 

143.00 

56.84 

1154.5 

6357.4 

101.47 

130.58 

11.38 

47.93 

71.38 

73.09 

145.00 

58.84 

1333.0 

6335.9 

97.43 

138.78 

13.66 

49.30 

l70. 10 

73.43 

147.00 

60.84 

1333.3 

6436.3 

99.53 

138.94 

13.30 

49.84 

69 . 46 

71.68 

149.00 

63.84 

1435.4 

6538.4 

103.37 

139.33 

13.94 

50.48 

^8.83 

70.83 

151.00 

64.84 

1518.3 

6631 .3 

105.16 

139.88 

14.44 

50.98 

^8.33 

70.05 

153.00 

66.84 

1601.8 

6704.8 

108.35 

130.60 

14.80 

51.34 

37 . 96 

69.35 

1 55 . 00 

68.84 

1760.7 

6863.7 

116.17 

133.76 

15.16 

51.70 

b7.60 

68.03 

157.00 

70.84 

1908.8 

7011.8 

130.83 

133.97 

15.80 

53.34 

36.96 

66.78 

159.00 

73.84 

1970.6 

7073.5 

119.89 

133.51 

16.44 

53.98 

36 . 33 

66 . 36 

160.40 

74.34 

3001.5 

7104.4 

119.44 

133.39 

16.76 

53.30 

;S6 . 00 

66.00 


*EOR 


\ 



Figure 7 

Samole Problem CALCOMP Plots 






Figiirp 7 (continued) 
Sample Problem CALfOMP Plots 
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Fiaure 7 (concluded) 
Sample Problem CALCOMP Plots 








APPENDIX A 


FORTRAN PROGRAM LISTING 


A complete FORTRAN source program listing is presented in the following 
pages. It starts with the main routine (UPSTAG) and is followed by the 
subroutines arranged in alphabetical order. There are a total of 481 cards in 
UPSTAG. The total program including subroutines contains 805 cards. 


- A-1 



uj nj 


1 

I 

ro 

I 


»DECK UPSTAG ( thpuT OUTPUT, TAPES - INPUT , TAPES - OUTPUT ) 

r nSSTAr^Tq^rPROGRAN FOR POST-FLIGHT ANALYSIS OF THE PITCH. YAU 

C and ROLL DISTURBANCES AND CONTROL FUEL CONSUMPTION FOR AN ON-OFF 

? DrSpirnN CONTROLLED UPPER STAOE ROCKET UEHICLE. 

® SnSnON^PUARIEOe 9),VUAR<2»e.9>.RVAR(B0O,9),TUAR(a08,9) 

ATMENSION nTITLE(8)»LTITL(16)»P(200)>T(200)»U(200),DPR(200)» 

DIMENSION NTllLtlB^^L 1 . 

?T?200) UP(200 KXCG(200),fiIXX(200),0IVY(2M^^^ 

Q ( 300 ). ALPH( 300), BETA! 300 ),CUAR( 200 > 10 ), BRACK ( 200 ) 

READI^sIsao! I0PT1,I0PT2,I0PT3,I0PTA,I0PT5,I0PTG 

C READ TUO CARDS OF ARBITRARY LABELING 

READ( 5.860 ) ( LTI TL ( J ) . J- 1 , 16 ) 

C READ TIPIE PARAMETERS 

READ( 5.850) TO, TBO. TSTEP. TCOAST. TSTEPC 
C INITIALIZE ALL COUNTERS. INDICES, AND SUMS 

AINT«0.0 
BINT-0.0 
CINT-0.0 
KKK-0 
LLL-0 
IC-0 


IT = 0 

TP-TO 

ITO-TO + 1 . 

ITCO-TBO+1 . 

MPHD-1 

MTT-1 

MQ-1 

MA-1 

MB-1 

MDP-1 

MDY-1 

MDR-1 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 



HR-1 36 

MP-1 37 

MT-1 38 

nu-1 39 

nUP-1 48 

nix -1 41 

HIY-1 43 

f1CG-l 43 

HTD-1 44 

HPSD*1 45 

IF (I0PT5 .LT. 0 ) GO TO 90 46 

C READ IH CONSTANTS 47 

READ( 5,850) XT, XC, ZC , RC, UFUELI , CNAS, XCP, ETA, CN3, DCP 48 

ETA«ETA/57.3 49 

C READ AERODYNAMIC PARAMETER TABLES Q, ALPHA, BETA 50 

READ( 5,880) NQ, ( Q (I ) , I > 1 , NO ) 51 

READ( 5,880 ) NALPH, ( ALPH( I ), I»1,NALPH ) 53 

READ( 5,880) NBETA, ( BETA ( I ) , I » 1 , NBETA ) 53 

C READ IN BOOSTER THRUST AND PROPELLANT REMAINING TIME HISTORIES 54 

READ( 5,840 ) NTT , ( TT ( I ) , I « 1 , NTT ) 55 

READ( 5,840 ) NNUP , ( UP ( I ) , I - 1 , NNUP ) 56 

C CONUERT PROPELLANT REMAINING TO PER-CENT UEIGHT REMAINING 57 

UTIGN»UP(3) 58 

DO 10 I«3,NNUP,3 59 

10 UP(I)«(UTIGN-UP(I))*100.0/UTIGN 60 

C READ IN MASS PROPERTIES UERSUS PERCENT BOOSTER PROPELLANT CONSUMED 61 

READ( 5,840 ) NUX, ( AIXX ( I ) , I - 1 , NUX ) 63 

READ( 5,840 ) NUY , ( AI YY ( I ), I - 1 , NUY ) 63 

READC 5,840 ) NXCG, ( XCG ( I ) , I - 1 , NXCG ) 64 

C TEST FOR INPUT OF ANGULAR ACCELERATION UALUES 65 

IF (lOPTl .GT. 0 ) GO TO 30 66 

C READ IN ANGULAR ACCELERATIONS 67 

READ( 5,840) NTHE , ( THEDD ( J ) , J « 1 , NTHE ) 68 

READ( 5,840 ) NPSI, (PSIDD( J ), J-1,NPSI ) 69 

READ( 5,840) NPHI, (PHIDD( J ), J-1,NPHI ) 70 



GO TO 90 


READ RATE TRACE SLOPE FACTORS 

READC 5,850) XKTHE , XKPSI , XKPHI , XfCTP , XKTY , XKTR 
READ IN SLOPE ANGLES 

READ( 5,840) NTHE, ( THEDD( I ), I » 1 ,NTHE ) 

READ( 5!840 ) NPS I , ( PSI DD ( I ) , I • 1 , NPS I ) 

READ( 5,840 ) NPHI , ( PH I DD ( I ) , I “ 1 , NPHI ) 

DO 30 I-a,NTHE,E 

THEDD( I )»THEDD( I )/57.3 

DO 40 I*2,NPSI,2 

PSIDD( I )»PSIDD( I )/57.3 

DO 50 I“2,NPHI,2 

PHIDDd )=PHIDD(I )/57.3 

CALCULATE ANGULAR ACCELERATIONS 

DO 60 I»2,NTHE,2 

THEDD( I ) = XKTHE*XKTP*TAN(THEDD( I ) ) 

DO 70 I=2,NPSI,2 

PSIDD( I ) = XKPSI*XKTY*TAN(PSIDD( I ) ) 

DO 80 I«2,NPHI,2 

PHIDDd ) = XKPHI*XKTRTTAN(PHIDDd ) ) 

TEST FOR INPUT OF ROLL IMPULSE TABLE 
READ( 5,840 ) NR I , ( R I (I ) , I * 1 , NR I ) 

TEST FOR INPUT OF DELTA RATES 
IF (I0PT3 .LE. 0 ) GO TO 130 
READ( 5,840 ) NDPR , ( DPR (I ) , I d , NDPR ) 

READ( 5,840 ) NDVR , ( DVR (I ) , I « 1 , NDYR ) 

READ( 5,840) NDRR , ( DRR (I ), I d , NDRR ) 

DPRATE=0.0 
DYRATE*0.0 
DRRATE-0.0 
DO 100 I*2,NDPR,2 
DPRd )»DPRATE + DPR(I ) 

DPRATEdPRd ) 

DO 110 1-2, NDYR, 2 


71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 



DYR( I)-DVRATE+DVR( I) 106 

110 DVRATE-DYR(I) 107 

DO 120 I-2,NDRR,2 108 

DRR( I )-DRRATE-^DRR( I ) 109 

120 DRRATE-DRR(I) 110 

C TEST FOR FUEL UEIGHT INPUT/CALCULATION 111 

130 IF (I0PT2 .GT. 0 ) GO TO 140 112 

C READ IN FUEL UEIGHT DATA 113 

READ( 5,840 ) NU, ( U ( I ) , I » 1 , NU ) 114 

GO TO 160 115 

C 116 

C FUEL UEIGHT CALCULATION 117 

140 READ( 5,850) DENS, COMFAC,PREREG, DUOL, V0LN2 118 

READ( 5,840) NNP, ( P ( I ) , I » 1 , NNP ) 119 

READ( 5,840 ) NNT, ( T( I ) , I = 1 ,NNT ) 120 

DO 150 I»2,NNT,2 121 

150 T( I )»T(I )+460.0 122 

C0NSTl»DENS»C0MFACl:U0LN2/PREREG 123 

C0NST2»DENS1:DU0L 124 

CALL C (P(2),T(2),CI ) 125 

IF (I0PT5 .GE. 0 ) GO TO 160 126 

READ( 5,850) RC,UFUELI , YHA,RIXX, YIYY,UFUELn,UFUELC,B0II1P,RIC 127 

GO TO 440 128 

C 129 

160 CONTINUE 130 

TTIHE*TP-TO 131 

C GET FUEL FUEL REMAINING ARRAY 133 

IF (I0PT2 .EQ. 0 ) GO TO 170 134 

CALL TBLU (NNP, PRES, TP, P, PIP ) 135 

CALL TBLU (NNT,TEP1P,TP,T,P1T ) 136 

CALL C (PRES,TEP1P CF ) 137 

DU-CONSTUC (P(2)*fEPiP )/(CI»T(2) )-PRES/CF ) + C0NST2»( (TEPIP/T(2 ) )-l .0 ) 138 

TUAR(KKK,9)-UFUELI-DU 139 

GO TO 180 140 



c 


170 


180 


C 

C 


I 

:p> 

I 

o^ 

I 


c 

c 


LOOKUP ROUTINE 


CALL TBLU ( NU , UUT , TP , U , NU ) 

I«?S5N'^tix’ivY’^CP IN TERNS OF TINE 

CfiLL^TBLU^(NNUP,UUP,TTinE.UP,nW 

mu !!lu?:Soiv^/utlp:Am:5iv) 

CfiU TBLU <NXCG,XXCGj^UUP,XCG^^ 

S^folSSSE^Nrut 

^Blu '(!IJs^;pJs“"d'/tp;JsidS:E1Jsd) 

CALL TBLU (NPHI ,PPHIDD,TP,PHIDD,«PHD ) 

COMPUTE THE SUM OF MOMENTS 

PUAR(KKK.7) = AAIYY1^TTHEDD/57.3 
VUAR(KKK 7)-AAIYY*PPSIDD/57.3 
R0AR(KKK,5 )=AAlXX*PPHIDD/57.3 

DEFINE XI AND X2 

COMPUTE^THRUST^MISAU COMPONENTS 

CALL TBLU (NTT,TTT,TTIME,TT,MTT ) 

TAIL TBLU ( NQ , QQ » TP , Q# MQ ) 

CALL TBLU ( NALPH , AALPH, TP, ALPH> MA ) 

CALL TBLU ( NBETA, BBETA, TP , BETA,MB ) 

IF (AALPH.NE.0. ) GO TO 190 
IF (BBETA.NE.0 . ) GO TO 190 
AERO»0. 

AB = 0 . 

PUAR(KKK,S)-0. 

YUAR(KKK,5)»0. 

GO TO 200 

COMPUTE TOTAL AERODYNAMIC ANGLE AND AERODYNAMIC MOMENTS 
190 TNA-TAN ( AALPH/57 . 3 ) 

TNB“TAN( BBETA/57 . 3 ) 


HI 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 



TANETA-SQRT(TNA*TNA + THB*TNB ) 176 

AB»57«3*ATAN(TANETA ) 177 

AERO- (CNAS*AB + CN3*AB**3 )*QO»(XXCG-XCP-DCP*ABS(AB ) )/12. 178 

PVAR(KKK,5)-AER0*TNA/TANETA 179 

VUAR(KKK,5 )«-AERO»TNB/TANETA 180 

200 PUAR(KKK,6)»PUAR(KIOC,7)-PUAR(KKK,5) 181 

VUAR(KKK,6)-VUAR(KKK,7)-YUAR(KKK,5) 182 

IF (TTT .GT. 0.) GO TO 210 183 

ETP-0.0 184 

ETV-0.0 185 

GO TO 220 186 

C 187 

210 ETP=(PUAR(KKK,7)-PUAR(KKK,5))/(TTT»X2) 188 

ETY=(YUAR(KKK,7)-YUAR(KKK,5) )/(TTT»X2 ) 189 

C COnPUTE TOTAL THRUST HISALIGNIIENT 190 

220 ET0T=SQRT(ETP**2+ETY**2) 191 

PUAR(K<K,8)-ETP*57.3 192 

YUAR(KKK,8)=ETY*57.3 193 

TUAR(KKK,4 )=ET0T*57.3 194 

IF (ETP) 230,270,310 195 

230 IF (ETY) 240,250,260 196 

240 ETOR-ATAN(ETY/ETP )*57. 3+180.0 197 

GO TO 350 198 

C 1 99 

250 ETOR-180.0 200 

GO TO 350 201 

C 202 

260 ETOR-ATAN(ETY/ETP )*57. 3+180.0 203 

GO TO 350 204 

C 205 

270 IF (ETY) 280,290,300 206 

280 ETOR-270.0 307 

GO TO 350 208 

C 209 

210 


290 ETOR-0.0 



GO TO 350 
C 

300 ETOR«90.0 
GO TO 350 
C 

310 IF (ETV) 330 , 330,340 

330 ETOR»ATAN(ETV/ETP )«57.3 + 360.0 
GO TO 350 
C 

330 ETOR=«0.0 
GO TO 350 
C 

340 ETOR = ATAN(ETY/ETP )*57.3 

350 TOAR(KKK,5)=ETOR 
C GET TIME ARRAY 

PUAR(KKK, 1 ) = TP 
YUAR(KKK, 1 )=TP 
ROAR(KKK, 1 )=TP 
TUAR(KKK, i )=TP 
PUAR(KKK,3)=TTIP1E 
YUAR(KKK,a )»TTIME 
RUAR(KKK,a )=TTIt1E 
TUAR(KKK,a)«TTIHE 
C GET BRACK ARRAY 

BRACKCKKK )» (XC-XXCG )*COS ( ETA )/ 13 . + ( ZC*SIN ( ETA ) )/13. 
C COMPUTE HYDROGEN PEROXIDE FUEL USED 

TUAR(KKK,8)*TUAR(1,9)-TUAR(KKK,9) 

C INTEGRATE MOMENTS IN PITCH AND YAU 

IF (KKK-1 .LE. 0 ) GO TO 370 

IF (I0PT6 .LE. 0 ) GO TO 360 

CALL TBLU (NDPR,DP,TP,DPR,MDP ) 

CALL TBLU ( NDYR, DY , TP, DYR ,MDY ) 

CALL TBLU (NDRR,DR,TP,DRR,MDR ) 
AINT»AAIYY*DP/(57.3*BRACK(KKK) ) 
BINT»AAIYY*DY/(57.3*BRACK(KKK) ) 


an 

aia 

ai3 

ai4 

ai5 

ai6 

ai7 

ai8 

319 

830 

331 

333 

333 

334 

335 

336 

337 

338 

339 

330 

331 
333 

333 

334 

335 

336 

337 

338 

339 

340 

341 
343 

343 

344 

345 
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CALL TBLU ( NRI , C INT , TP, RI , MR ) 346 

CINT-CINT+AAIXX«DR/(57.3*(RC/ia.0 ) ) 347 

GO TO 370 248 

C 249 

360 AINT-AINT+ABS( (PUAR(ICICK,7) + PUAR(ICKK-1,7) )«(TUAR(KKK,1 )-TUAR(KlCK-l, 350 
1 1 ) ) / ( a ♦ 0< BRACK ( KKK ) ) ) 351 

BINT-BINT+ABS( (VUAR(KKK,7)+YUAR(KKK-1,7) )»(TUAR(KKK,n-TUAR(KKK-l, 353 
11 ) )/(3.0»BRACK(KKK) ) ) 353 

CALL TBLU ( NRI , C INT , TP ,RI , MR ) 354 

CINT»CINT + ABS(6.0*(RUAR(KKK,5) + RUAR(KKK-1,5 ) )*(TUAR(KKK, 1 )- 355 

1 -TUAR(KKK-1,1) )/RC) 356 

C FIND TOTAL IMPULSE 357 

370 TUAR(KKK,7 )=AINT+BINT+CINT 358 

C COMPUTE THE SPECIFIC IMPULSE 359 

IF (TUAR(KKK,8) .NE. 0 ) GO TO 380 360 

TUAR(KKK,6)»0.0 361 

GO TO 390 363 

C 263 

380 TUAR(KKK,6)=TUAR(KKK,7 )/TUAR(KKK,8) 264 

390 CONTINUE 265 

PUAR(KKK,3 )=TTHEDD 266 

PUAR(KKK,4)=AAIYV 367 

PUAR(KKK,9)*AINT 268 

YUAR(KKK,3)=PPSIDD 269 

YUAR(KKK,4)=AAIYY 270 

YUAR(KKK,9)=BINT 271 

RUAR(KKK,3)=PPHIDD 272 

RUAR(KKK,4 )=AAIXX 273 

RUAR(KKK,6)»CINT 274 

RUAR(KKK,7)=XXCG 275 

TUAR(KKK,3)=TTT 276 

IF (IT .NE. 0 ) GO TO 400 277 

TP=ITO 378 

IT»1 379 

GO TO 160 380 
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400 IF (TBO-(TP+TSTEP ) .LE. 0 ) GO TO 410 
TP-TP+TSTEP 
GO TO 160 
C 

410 IF (TBO-TP .LE. 0 ) GO TO 430 
TP-TBO 
GO TO 160 
C 

420 AUGISP»TUAR(KKK,6 ) 

DO 430 1=1, KKK 

ROAR( I,8)=TUAR( I,7)/AUGISP 

430 RUAR(I,9)»TUAR(1,9)-RUAR(I,8) 

IF (I0PT5 .LE. 0 ) GO TO 580 
B0If1P»TUAR(KKK,7) 

UFUELt1»TUAR(KKK,9 ) 

UFUELC = RUAR(KKK,9 ) 

RIC=RUAR(KKK,6) 

YHA = BRACKOCKK ) 

RIXX=AAIXX 

VIYY=AAIVV 

440 CONTINUE 

TTIf1E = TP-TO 

LLL=LLL+1 

CUAR(LLL,1 ) = TP 

CUAR(LLL,2 ) = TTI(1E 

IF (I0PT2 ,GT. 0 ) GO TO 450 

CALL TBLU ( NU, UUT , TP , U, MU ) 

CUAR(LLL,9)»UUT 
GO TO 460 
C 

450 CALL TBLU ( NNP , PRES , TP . P, HP ) 

CALL TBLU (NNT,TEMP,TP,T,nT ) 

CALL C (PRES TEMP CF ) 

DU=CONSTl»( (P(2)»fEMP)/(CI*T(2) )-PRES/CF )+C0NST2«( (TEHP/T(2) 


281 

282 

283 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 

296 

297 

298 

299 

300 

301 

302 

303 

304 

305 

306 

307 

308 

309 

310 

311 

312 

313 

314 
0) 315 
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CUAR(LLL,9)»UFUELI-DU 31B 

460 CUAR(LLL,8)»UFUELI-CUAR(LLL>9) 317 

CUAR(LLL,7 )-UFUELH-CUAR(LLL,9 ) 318 

IF (I0PT3 .LE. 0 ) GO TO 470 319 

CALL TBLU ( NDPR, DP > TP, DPR , MDP ) 380 

CALL TBLU ( NDVR , DY, TP , DVR , HDV ) 381 

CALL TBLU ( NDRR , DR , TP, DRR , MDR ) 388 

PinP-YIYY*DP/(57.3»VnA) 383 

YinP = YIYY*DY/(57.3*YriA) 384 

RinP>'RIXX*DR/(57.3l:(RC/ia.0) ) 385 

C0AR(LLL,3) = PIHP + YinP + RII1P 386 

CUAR(LLL,4)=CUAR(LLL,3)+B0IHP 387 

GO TO 480 388 

C 389 

470 CALL TBLU ( NRI , T IMP , TP, RI , HR ) 330 

CUAR(LLL,4 )=BOIHP+TinP-RIC 331 

CUAR(LLL,3)=CUAR(LLL,4 )-BOIHP 338 

480 IF (CUAR(LLL,7) ,NE. 0.) GO TO 490 333 

CUAR(LLL,5)=0.0 334 

GO TO 500 335 

C 336 

490 CUAR(LLL,5)=CUAR(LLL,3)/CUAR(LLL,7) 337 

500 IF (CUAR(LLL,8) .NE. 0.) GO TO 510 338 

CUAR(LLL,6 )"0.0 339 

GO TO 580 340 

C 341 

510 CUAR(LLL,6)=CUAR(LLL,4 )/CUAR(LLL,8) 348 

580 IF (IC .GT. 0 ) GO TO 530 343 

TP=ITCO 344 

IC=1 345 

GO TO 440 346 

C 347 

530 IF (TCOAST-(TP+TSTEPC ) .LE. 0.) GO TO 540 348 

TP»TP+TSTEPC 349 

GO TO 440 350 
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540 IF (TCOAST-TP .LE. 0.) GO TO 550 
TP^TCOAST 
GO TO 440 

550 AUGISP»CUAR(LLL,5 ) 

DO 570 I-1,LLL 

IF (AUGISP .NE. 0. ) GO TO 560 

CFUELC-0.0 

GO TO 570 

560 CFUELC»CUAR( 1,3 )/AUGISP 
570 CUAR(I,10)»LIFUELn-CFUELC 

IF (I0PT5 .LT. 0 ) GO TO 630 
580 XKKK=KKK 
XP=XKKK/50. 

NP=KKK/50 

XNP=NP 

IF(XP.GT.XNP)NP=NP+1 
NLL = 0 
NPP = 1 

590 IF (NP-NPP .GT. 0 ) GO TO 600 
NFL*NLL+1 
NLL=KKK 
GO TO 610 

600 NFL=NLL+1 
NLL»NLL+50 

610 URITEC 6,690) NPP 

URITE( 6,870) (LTITU J), J-1,16) 

URITE( 6,700) 

URITE( 6,710) ((PUAR(I,J),J«1,9),I»NFL,NLL) 

URITE( 6,730) NPP 

URITE( 6,870) (LTITL( J ), J»l, 16) 

URITE( 6,730) 


351 

353 

353 

354 

355 

356 

357 

358 

359 

360 

361 
363 

363 

364 

365 

366 

367 

368 

369 

370 

371 
373 

373 

374 

375 

376 

377 

378 

379 

380 

381 
383 

383 

384 

385 
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URITE( 6,710) ((YUAR(I,J),J-1,9),I“NFL,NLL) 386 

URITE( 6,740) NPP 387 

URITE( 6,870) (LTITL(J),J«1,16) 388 

URITE( 6,750) 389 

URITE( 6,760) ((RUARd, J), J-1,9),I“NFL,NLL) 390 

URITE( 6,770) NPP 391 

URITE( 6,870) (LTITL(J),J»1,16) 398 

MDTTPf A 1 '^Q'5 

URITE( 6'880) ((TUAR(I,J),J»1,9),I»NFL,NLL) 394 

IF (NP-NPP .LE. 0 ) GO TO 680 395 

NPP»NPP+1 396 

GO TO 590 397 

C 398 

680 CONTINUE 399 

IF (I0PT5 .LE. 0 ) GO TO 670 400 

630 XLLL»LLL 401 

XC=XLLL/50. 408 

NC«LLL/50 403 

XNC=NC 404 

IF(XC.GT.XNC)NC=NC+1 405 

NCLL=0 406 

NCPP=1 407 

640 IF (NC-NCPP .GT. 0 ) GO TO 650 408 

NCFL=NCLL+1 409 

NCLL=LLL 410 

GO TO 660 411 

C 418 

650 NCFL»NCLL+1 413 

NCLL-NCLL+50 414 

660 URITE( 6,790) NCPP 415 

URITE( 6,870) (LTITU J),J«1,16) 416 

URITE( 6,800) 417 

URITE( 6,810) ((CUAR(I,J),J»1,10),I«NCFL,NCLL) 418 

IF (NC-NCPP .LE. 0 ) GO TO 670 419 

NCPP-NCPP+1 4B0 
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GO TO 640 4ai 

c 4aa 

670 CONTINUE 4a3 

IF (IOPT4 .LE. 0 ) GO TO 680 ABA 

READ( 5,860 ) ( NT I TLE ( J ) , J = 1 , 8 ) 485 

NQPL-0 486 

IF(Q(8).GT.0. )NQPL»1 487 

READ( 5,840 ) IROLL , SFT , SFPVd , SFET , SFRN 488 

CALL CURUE (KKK,NQPL, IROLL, NTITLE,5FPVf1, SFRN, SFET, SFT ) 489 

680 CONTINUE 430 

STOP 431 

C THIS ENDS THE UPSTAG EXECUTION. ONLY FORMAT STATEMENTS REMAIN. 438 

690 FORMAT ( IHl , 85X, 86HUPPER STAGE MOMENT ROUTINE, 14X,6HPAGE »,I8,8H A 433 
1/85X,31H PITCH UARIABLES TIME HISTORIES) 434 

700 FORMAT (//4X,6HFLIGHT.8X,5HSTAGE,aX,7HANGULAR,aX,6HM0MENT,4X,4HAER 435 
10,3X,8HMISALIGN,3X,5HT0TAL,3X,6HTHRUST/8(4X,4HTIME ),4X,5HACCEL,aX, 436 
88H0F INER. ,8X,3(6HM0MENT,3X ),6HMISALN, 1X,7HIMPULSE/5X,3HSEC,5X,3HS 437 
3EC,3X,7HDEG/S/S, 1X,8HSLUG-FT8,3X,3(5HFT-LB,4X),3HDEG,4X,6HLB-SEC// 438 
4 ) 439 

710 FORMAT (3X,F7.a,lX,F6.a,2X,F6.a,8X,F8.1,lX,F8.a,lX,F8.a,lX,F8.a,lX 440 

1,F7.4,1X,F7.8) 441 

780 FORMAT ( IHl , 85X, 86HUPPER STAGE MOMENT ROUTINE, 1 4X, 6HPAGE «,I8,8H B 448 
1/85X,31H VAU UARIABLES TIME HISTORIES) 443 

730 FORMAT ( //4X, 6HFLIGHT , 8X , 5HSTAGE , 8X, 7HANGULAR, 8X, 6HM0MENT, 4X , 4HAER 444 
10,3X,8HMISALIGN,3X,5HT0TAL,3X,6HTHRUST/8(4X,4HTIME ),4X,5HACCEL,aX, 445 
88H0F INER. ,SX,3(6HM0MENT,3X ),6HMISALN,1X,7HIMPULSE/5X,3HSEC,5X,3HS 446 
3EC,3X,7HDEG/S/S, IX, 8HSLUG-FT8, 3X, 3(5HFT-LB, 4X ), 3HDEG, 4X, 6HLB-SEC// 447 
4 ) 448 

740 FORMAT ( IH 1 , 85X , 86HUPPER STAGE MOMENT ROUTINE, 14X,6HPAGE -,I8,8H C 449 
1/85X,31H ROLL UARIABLES TIME HISTORIES) 450 

750 FORMAT ( //4X, 6HFLIGHT , 8X, 5HSTAGE , 8X, 7HANGULAR, 1 X , 6HM0MENT, 4X , 5HT0T 451 
lAL, 13X,4HC.G. ,4X, 18H CALCULATED FUEL /8(4X, 4HTIME ),4X,5HACCEL,8X, 458 
87H0F INER,3X,6HM0MENT,8X,7HIMPULSE,3X,5HP0INT,3X,7H CONSUM, 3X , 7HRE 453 
3MAIN /5X,3HSEC,5X,3HSEC,3X,7HDEG/S/S, 1 X , 8HSLUG-FT8 , 8X , 5HFT-LB , 3X, 6 454 

4HLB-SEC,4X,4HINCH,7X,8HLB,7X,8HLB//) 455 


I 
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760 FORMAT (3X#F7»2*1X»F6»2>3X»F6»2»2X>F8.1>1X#F8»E# lX«F7*2«lXfF8*3«lX 
1 F8 3 IX F8 3) 

770 FORMAT ( IHl ) 25X, 26HUPPER STAGE MOMENT ROUTINE, 14X,6HPAGE -,I2,2H D 
1/25X,31H TIME HISTORIES OF SYSTEM ) 

780 FORMAT ( //4X, 6HFLIGHT , 2X, 5HSTAGE , 1 IX, 5HT0TAL, lOX, 8HSPECIFIC, 2X,5HT 
10TAL,2(4X,4HFUEL)/2(4X,4HTIME),3X,6HTHRUST,2X,7HMISALGN,2X,5HLAMDA 
2,2X,2(7HIMPULSE,2X),GHC0NSUM,2X,6HREMAIN/5X,3HSEC,5X,3HSEC,5X,2HLB 
3,6X,3HDEG,5X,3HDEG,5X,3HSEC,4X,6HLB-SEC,4X,2HLB,6X,2HLB//) 

790 FORMAT ( IHl , 2BX, 26HUPPER STAGE MOMENT ROUTINE, 14X,GHPAGE -,I2,2H E 
1/25X,31H COAST UARIABLES TIME HISTORIES) 

800 FORMAT ( //, 4X, 6HFLIGHT, 2X, 5HSTAGE, 3X, 5HC0AST, 3X, BHTOTAL, 2X, 5HC0AST 
1,,2X,5HT0TAL,3X,5HC0AST,3X,5HT0TAL,2X,14HFUEL REMAINING, /5X, 4HTIME 
2,3X,4HTIME,3X,2(7HIMPULSE,1X),1X,2(3HISP,4X),4HFUEL,4X,4HFUEL,3X,1 
34H MEAS CALC /5X,3HSEC,5X,3HSEC,3X,2(6HLB-SEC,2X), 1X,2(3HSEC,4X 
4),2HLB,7X,2HLB,6X,2HLB,6X,2HLB//) 

810 FORMAT (3X,F7,2,1X,F6.2,1X,F7.1,1X,F7.1,1X,F6.2,1X,F6.2,1X,F7.2,1X 
1 F7.2 IX F7.2 IX F7.E) 

820 FORMAT ( 3X, F7 ! 2 , 1 X , F6 . 2, 2X, F7 . 0, IX, F7 . 4 , IX, F7 . 2 , IX, F7 . 2 , IX, F8 . 2, IX 


1,F6.2,2X,F7,2) 


830 FORMAT 
840 FORMAT 
850 FORMAT 
860 FORMAT 
870 FORMAT 
880 FORMAT 
END 


(1015) 

(I5/(6E10.3) ) 
(6E10.3) 
(8A10) 
(1X,8A10) 
(I5/(7E10.3) ) 
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471 
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474 
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^DECK C 1 

CIIRDAIITTNP f- ( TAT ^'AT ^ S 

THIS SUBROUTINE COMPUTES THE COMPRESSIBILI TV FACTOR FOR DRV 3 

NITROGEN FOR 50-100 DEG F, AND 400-3000 PSIA 4 

PP-PAT-1400. 5 

CAT- .9977+1 .657738E-5*PP+1 .E64881E-8*PP»PP+4.5833E-4*(TAT-5aa. ) 6 

RETURN 

END 8 


CK CURUE 1 

SUBROUTINE CURUE (NP, NQ, IROLL,NTITLE, DV, DVR, DVET, DX ) 3 

THIS SUBROUTINE PERFORMS THE CALCOMP PLOT SETUP FOR PITCH, YAU 3 

AND ROLL MOMENTS AND THRUST MISALIGNMENT. 4 

DIMENSION NTITLE(8),X(400),V(400) 5 

COMMON PUAR(200,9),YUAR(200,9),RUAR(a00,9),TUAR(E00,9) 6 

DATA Zl,Z2,Z3/'0. 42,0.07,0. 14/ 7 

SET UP INITIAL AXIS UARIABLES 8 

9 

VF“-2.*DY 10 

NDX-DX 11 

CALL PLOTS (5HCAL19,0,4HPLOT) 12 

DO 10 J“1,NP 13 

10 X(J)«PUAR(J,1 ) 14 

C PLOT PITCH MOMENTS 15 

CALL SYMBOL (0 . 5 , 9 . 9, 0 . 10, NTITLE, 0 . , 80 ) 16 

CALL SYMBOL ( 0 . 5 , 9 . 5, 0 . 10, 45H PITCH AND YAU DISTURBING 17 

1 MOMENTS, 0. ,45 ) 18 

CALL SYMBOL(3.5,9.3,0.07,12HTOTAL MOMENT, 0 ., 1 2 ) 19 

CALL PL0T(3.4,9.3,3 ) 20 

CALL PL0T<2.4,9.3,2 ) 21 

CALL PL0T<2.4, 9. 15, 3 ) 22 

CALL PL0T(2.6, 9. 15,2 ) 23 

CALL PL0T(2.7,9. 15,3) 24 

CALL PL0T(2.9, 9. 15,2 ) 25 



CALL PLOTO. 0,9. 15,3 ) 86 

CALL PL0T(3.a,9.15,a) 87 

CALL PL0T(3.3,9. 15,3 ) 88 

CALL PL0T(3.4,9. 15,8 ) 89 

CALL SYMBOL(3.5,9.15,0.O7,26HTHRUST MISALIGNnENT flOMENT ,0.,26) 30 

CALL PLOT(B.4,9.0,3 ) 31 

CALL PLOT(a.8,9.0,8) 38 

CALL PLOT(8.90,9.0,3) 33 

CALL PLOT(8.95,9.0,a) 34 

CALL PLOTO.05,9.0,3) 35 

CALL PLOT(3.4,9.0,8) 36 

CALL SVnBOL(3.5,9.0,0.07,18HAERODYNAfniC MOriENT ,0.,18) 37 

NTI=X(1)/DX 38 

XF=NTI»NDX 39 

CALL PLOT (l.,7.,-3) 40 

CALL AXIS (0.,-S.,80HFLIGHT TIME SECONDS ,-80,6. ,0. ,XF,DX ) 41 

DO 20 J-1,NP 48 

X(J )“(X( J)-XF )/DX 43 

20 V( J )-PUAR(J,7) 44 

CALL AXIS (0.,-8.,19HPITCH MOMENT FT-LBS, 19,4. ,90. ,VF,DY ) 45 

CALL DASH ( X, Y, NP, Z1 , Z1 , 0 . , 1 . , DY,0, 7 . , 3. ) 46 

IF (NQ.EQ.0) GO TO 50 47 

DO 30 J-1,NP 48 

30 Y( J )»PUAR( J,6) 49 

CALL DASH ( X, Y, NP, Z3, Z3, Z8, 1 . , DY, 0, 7 . , 3. ) 50 

DO 40 J-1,NP 51 

40 Y( J )-PUAR( J,5 ) 52 

CALL DASH ( X, Y, NP, Z1 , Z2, Z3, 1 . , DY, 0, 7 . , 3 . ) 53 

C PLOT YAU MOMENTS 54 

50 CALL PLOT (0.,-5.,-3) 55 

CALL AXIS (0. ,-2. ,20HFLIGHT TIME SECONDS , -20, 6 . , 0 . , XF, DX ) 56 

CALL AXIS (0.,-2.,17HYAU MOMENT FT-LBS, 17, 4 . , 90. , YF, DY ) 57 

DO 60 J-1,NP 58 

60 Y( J )-YUAR( J,7) 59 

CALL DASH ( X, Y, NP, Z1 , Z1 , 0 . , 1 . , DY, 0, 7 . , 3 . ) 60 


IF (NQ.EQ.0) GO TO 90 61 

DO 70 J-1,NP 63 

70 V< J )-YUAR( J,6 ) 63 

CALL DASH ( X, Y, HP, Z3, 23 , 23, 1 . , DV, 0, 7 . , 3 . ) 64 

DO 80 J-1,NP 65 

80 Y( J )-YUAR( J,5 ) 66 

CALL DASH ( X, Y, NP, 21 , 23, 23, 1 . , DY, 0, 7 , , 3 . ) 67 

90 CALL PLOT (-l.,-3.,-3) 68 

CALL PLOT <13. ,0., -3) 

CALL PLOT (0.,0.,999) 70 

C CHECK AMD PLOT ROLL nOMENTS 71 

IF (IROLL.EQ.0) GO TO 110 73 

YFR=-4.*DYR 73 

CALL PLOTS (5HCAL19,0,4HPLOT) 74 

CALL SYMBOL ( 1 . , 9 . 8, 0 . 10, NTITLE, 0 . , 80 ) 75 

CALL SYMBOL (1 . , 9 . 5, 0 . 10, 35H ROLL DISTURBING MOMENT, 0. 76 

1,35) 77 

CALL PLOT <l.,5.,-3) 78 

CALL AXIS (0. ,-4. ,30HFLIGHT TIME SECONDS ,-30,6. ,0. ,XF,DX ) 79 

CALL AXIS (0.,-4.,39HROLL DISTURBING MOMENT FT-LBS, 39, 8. , 90. , YFR, D 80 
IVR) 81 

DO 100 J“1,NP 83 

100 Y( J )-RUAR< J,5 ) 83 

CALL DASH ( X, Y, NP, 21 , 21 , 0 . , 1 . , DYR,0, 7 . , 5 . ) 84 

CALL PLOT (11., -5., -3) 85 

CALL PLOT (0.,0.,999) 86 

C PLOT THRUST MISALIGNMENTS 87 

110 CALL PLOTS (5HCAL19,0,4HPLOT) 88 

YFET=-DYET 89 

CALL SYMBOL (1 . , 9 . 8, 0 . 10, NTITLE , 0 . , 80 ) 90 

CALL SYMBOL ( 1 . , 9 . 5 , 0 . 1 0, 33H THRUST MISALIGNMENT, 0 ., 33 91 

1 ) 93 

CALL PLOT (l.,7.,-3) 93 

CALL AXIS (0. ,-l . ,30HFLIGHT TIME SECONDS , -30 , 6 . , 0 . , XF, DX ) 94 

CALL AXIS (0. , -1 . , 19HPITCH COMPONENT DEG, 19, 3 . , 90 . , YFET, DYET ) 95 



DO 120 J-1,NP 96 

120 V<J)-P0AR(J,8) 97 

CALL DASH ( X, V, NP, Z1 , Z1 . 0 . , 1 . > DVET, 0, 7. , 3. ) 98 

CALL PLOT (0.,-3.,-3) 99 

CALL AXIS (0, .-l.,20HFLIGHT TIME SECONDS ,-20,6. ,0. ,XF,DX ) 100 

CALL AXIS (0. ,-l ♦ , 17HYAU COflPONENT DEG, 17, 2. , 90 . , YFET, DYET ) 101 

DO 130 J-1,NP 102 

130 Y( J )-YOAR( J,8) 103 

CALL DASH (X, Y,NP,21,Z1 ,0. , 1 , ,DYET,0,7. ,3. ) 104 

CALL PLOT (0,,-4.,-3) 105 

CALL AXIS (0.,0. ,20HFLIGHT TIME SECONDS ,-20,G. ,0. ,XF,DX ) 106 

CALL AXIS I0.,0.,9HTOTAL DEG, 9, 2. ,90. ,0. , DYET ) 107 

DO 140 J»1,NP 108 

140 Y(J)=TUAR(J,4) 109 

CALL DASH ( X, Y, NP, 21 , Z1 , 0 . , 1 . , DYET, 0, 7. , 3 . ) 110 

CALL PLOT (11. ,0., -3) 111 

' CALL PLOT (0.,0.,999) 112 

> RETURN 113 

END 114 
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ooo 


I 
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0 
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tDECK 


10 

E0 

C 

30 


40 


DASH 1 
SUBROUTINE DASH ( X , V , NP , Z1 , Z3, SPACE , XSCALE , YSCALE > LSVMB , XL IM , VL IN ) a 
SYMBOLS, DASHED, DASHED-DOT LINES OR SOLID LINES UITH OR UITHOUT 3 
SYMBOLS BASED ON A SET OF SEQUENTIAL POINTS GIUEN IN 4 
THE INPUT 'X' ABSCISSA ARRAY AND THE 'V' ORDINATE ARRAY 5 
DIMENSION X(1 ),Y(1 ) 6 
DO 10 I-l.NP 7 
XA“X( I )/XSCALE 8 
YA-Y( I )/VSCALE 9 
IF (ABS(XA).GT.XLIM) GO TO 10 10 
IF (ABS(YA).GT.YLIM ) GO TO 10 11 
CALL PLOT (XA,YA,3) 13 
GO TO 30 13 
CONTINUE H 
IF (SPACE) 330,310,30 15 
THIS SUBROUTINE PLOTS A CALCOMP PLOT UITH A UIDE UARIETY OF 16 
K = 0 17 
PIE“1.5708 18 
Z»Z1 19 
ZB=Za 20 
IF (za .GT. 0. ) GO TO 40 31 
ZB=Z1 23 
ZD=Z 23 
LZ»0 24 
SL=0. 25 
NF-NP-1 26 
DO 300 J=1,NF 27 
XA»X( J )/XSCALE 28 
IF (ABS(XA)-XLIM ,GT. 0.) GO TO 300 29 
XB-X( J + 1 )/XSCALE 30 
IF (ABS(XB)-XLIM .GT. 0.) GO TO 300 31 
YA-Y( J)/YSCALE 32 
IF (ABS(YA)-YLIM .GT. 0.) GO TO 300 33 
YB-Y( J + 1 )/YSCALE 34 
IF (ABS(YB)-YLin .GT. 0.) GO TO 300 35 



DV«YB-YA 36 

DX-XB-XA 37 

IF (DX .HE. 0. ) GO TO 80 38 

IF (DY) 50>60,70 39 

50 TH--Pia 40 

GO TO 90 41 

£0 TH-0. 42 

GO TO 90 43 

70 TH-PI2 44 

GO TO 90 45 

80 TH«ATAN(DY/DX) 46 

90 DX«XB-XA 47 

DY»YB-YA 48 

DZ»SQRT(DX*DX+DY*DY) 49 

C TEST TO SEE UHAT IS GOING ON 50 

IF (K) 100,180,220 51 

100 K»1 52 

c I 3 C PA TF F 1 

IF (DZ-SPACE) 110,120,150 54 

C SPACE IS LARGER THAN DZ 55 

110 SL-SL-DZ 56 

CALL PLOT (XB,YB,3) 57 

GO TO 300 58 

C NEXT POINT IS EXACTLY ONE SPACE 59 

120 K-0 60 

IF (LZ .NE. 0 ) GO TO 130 61 

ZD-ZB 62 

LZ-1 63 

GO TO 140 64 

130 ZD«Z 65 

LZ-0 66 

140 SL-0. 67 

CALL PLOT (XB,YB,3) 68 

GO TO 300 69 

C NEXT POINT MORE THAN ONE SPACE AUAY 70 



150 XA-XA+SPACEl:COS(TH) 71 

YA-YA+SPACE»SIN(TH) 72 

IF (ABS(XA)-XLItl .GE. 0.) GO TO 300 73 

IF (ABS(YA )-YLIf1 .GE. 0.) GO TO 300 7H 

K = 0 75 

IF (LZ .HE. 0 ) GO TO 160 76 

ZD-ZB 77 

LZ*1 78 

GO TO 170 79 

160 ZD'Z 80 

LZ»0 81 

170 SL=0. 82 

CALL PLOT (XA,YA,3) 83 

GO TO 90 8^ 

C K=0 LINE BEING DRAUN ZD LENGTH NOT DRAUN RESUME AS IS LINE STARTING 85 

180 IF (DZ-ZD) 190,200,210 86 

C LINE GOES AT LEAST TO NEXT POINT 87 

190 K»0 88 

ZD»ZD-DZ 89 

CALL PLOT (XB,YB,2) 90 

GO TO 300 91 

C LINE ENDS AT NEXT POINT 92 

200 K=-l 93 

SL>SPACE 94 

ZD=0. 95 

CALL PLOT (XB,YB,2) 96 

GO TO 300 97 

C LINE ENDS BEFORE NEXT POINT 98 

210 K-1 99 

SL-SPACE 100 

XA»XA+ZD*COS(TH) 101 

YA-YA+ZD*SIN(TH) 102 

IF (ABS(XA)-XLIM .GE. 0.) GO TO 300 103 

IF (ABS(YA)-YLin .GE. 0.) GO TO 300 104 

CALL PLOT (XA,YA,2) 105 



GO TO 90 107 

C K«1 IS IN SPACE 108 

H20 2D-0. 109 

IF (DZ-SL) 230,240,270 110 

230 K-1 111 

SL-SL-DZ iia 

CALL PLOT (XB,YB,3) 113 

GO TO 300 IM 

C SL-DZ 115 

240 K=0 116 

IF (LZ .NE. 0 ) GO TO 250 117 

ZD=ZB 118 

LZ=1 119 

GO TO 260 120 

250 ZD=Z 121 

LZ=0 122 

260 CALL PLOT (XB,YB,3) 123 

GO TO 300 124 

C SL IS LESS THAN DZ 125 

270 K=0 126 

IF (LZ .NE. 0 ) GO TO 280 127 

ZD»ZB 128 

LZ»1 129 

GO TO 290 130 

280 ZD=2 131 

LZ=0 132 

290 XA = XA + SLl:COS(TH ) 133 

YA«YA+SL*SIN(TH) 134 

IF (ABS(XA)-XLIH .GE. 0.) GO TO 300 135 

IF (ABS(YA)-YLin .GE. 0.) GO TO 300 136 

SL»0. 137 

CALL PLOT (XA,YA,3) 138 

GO TO 90 139 

300 CONTINUE 140 



c 


c 


GO TO 370 

STRAIGHT LINE PLOT OPTION 
310 DO 330 J-I>NP 
XA-X(J )/XSCALE 
VA-Y( J)/VSCALE 

IF (ABS(XA)-XLin .GT. 0.) GO TO 330 
IF (ABS(YA)-YLIH ,GT» 0.) GO TO 330 
CALL PLOT (XA,YA,3) 

330 CONTINUE 
GO TO 370 

PLOT SYMBOLS ON LINE NO LINE IF LYSMB IS NEGATIUE 
330 NSM=IABS(LSYMB ) 

IF (LSYMB .LT. 0 ) GO TO 340 

K s -2 

GO TO 350 


340 K*-l 

350 DO 360 J=1,NP 
XA=X( J )/XSCALE 
VA = Vf M/V^TAIF 

IF (ABS(XA)-XLIM .GT. 0.) GO TO 360 
IF (ABS(YA)-YLIM .GT. 0.) GO TO 360 
CALL SYMBOL ( XA , YA , 0 . 07 , NSM, 0 . 0, K ) 
360 CONTINUE 
370 CALL PLOT (0.,0.,3) 

RETURN 

END 
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oooooooo 


tDECK TBLU 

SUBROUTINE TBLU (NT,Y,X,T,H) 

SINGLE TABLE LOOKUP SUBROUTINE 
NT • NUMBER OF UALUES IN ARRAY 
Y - RETURNED ORDINATE 
X - ABSCISSA UALUE CALLED 

T - INPUT TABLE OF ALTERNATING ABSCISSAS AND ORDINATES 
ORDINATES MUST BE HONOTONICALLY INCREASING 
M • PREUIOUS INDEX USED IN THIS TABLE LOOKUP 
THIS INDEX GETS CHANGED TO CURRENT UALUE 
DIMENSION T( 1 ) 

10 IF (T(M)-X) 50,ae>30 
a0 Y»T(M+1 ) 

RETURN 

30 IF (T(l )-X.LT.0. ) GO TO 40 
M-1 

GO TO ao 
40 M = f1-a 
GO TO 10 
50 MM = I1 + a 

IF (rin-NT-l.LE.0) GO TO 60 
n=NT-l 
GO TO a0 

60 IF (T(MM)-X.GT.0. ) GO TO 70 
fl'lin 

GO TO 50 
70 M»nn-a 

DT=T(MM)-T(M) 

IF (DT.NE.0. ) GO TO 80 
Y-T(M+1 ) 

RETURN 

80 DY«T(MH + 1)-T(M+1 ) 

DDT-X-T(H) 

Y»T(M + 1 )+DY*DDT/DT 
RETURN 
END 
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